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SUMMARY 


In  this  Report  the  closed-loop  control  system  for  the  ballistic  test 
vehicle  Black  Knight  is  described  and  an  account  is  given  of  the  theory  under¬ 
lying  its  design.  Attitude  control  of  the  vehicle  is  achieved  by  swivelling 
the  four  rocket  motor  combustion  chambers  in  response  to  gyroscope  and  ccmnand 
guidance  signals.  The  effects  of  propellent  sloshing,  aerodynamic  instability 
euid  flexibility  of  the  vehicle  stricture  are  among  the  factors  taken  into 
account  in  determining  transfer  fxaiotions,  relating  combustion  chamber  move¬ 
ments  to  gyroscope  signals,  viiich  satisfy  the  conditions  for  stable  fli^t. 

In  the  final  section  a  brief  account  is  given  of  ground  tests,  vfcich 
resulted  in  some  alterations  to  the  roll  stabilisation  system,  and  of  early 
flight  trial  results.  The  control  system  has  performed  very  satisfactorily 
in  all  flints  to  date. 
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1  INTRCiDaCTIQN 

•j 

Blaok  Kni^t  is  a  balliatio  test  vehicle  intended  for  research  into 
problems  associated  with  ballistic  missiles,  especially  that  of  hi^  speed 
re-entry  into  the  atmosphere.  The  develppnent  of  Black  Knight  is  also 
e^qjected  to  yield  usefiil  experience  in  the  design  and  operation  of  the  type 
of  control  system  used  in  ballistic  missiles. 

The  layout  of  Black  Khigjit  is  shown  in  Pig*1,  Pour  swivelling  rocket 
motors,  burning  kerosene  and  hydrogen  peroxide  (H.T.P,),  are  used  for  pro- 
pvilsion  and  attitude  control  of  the  vehicle.  Each  motor  has  a  sea-level 
thrust  of  4000  lb  and  can  be  rotated  througjh  HO®  about  a  single  axis  by  an 
electro-hydraulic  position  control  servo-mechanism.  Opposite  pairs  are 
deflected  together  to  provide  pitch  and  yaw  attitude  control  and  all  four 
differentially  for  roll  stabilisation.  The  motor  pivot  axes  are  so  inclined 
that  the  direction  of  thrust  makes  an  an^e  of  with  the  vehicle  centre 
line,  in  order  to  reduce  disturbances  ydiloh  migjit  arise  if  the  motor  thrusts 
were  xinequal,  partictilarly  near  all-burnt.  The  swivelling  combustion 
chambers  are  supported  by  a  oonposite  sti*uoture  which  transmits  thrust  and 
other  loads  to  the  base  of  the  H,T,P,  tank.  Propellents  are  pumped  in 
throu^  the  trxjnnion  mountings,  kerosene  throu^  the  outer  bearing  and  H,T.P, 
throu^  the  inner.  The  3  ft  diameter  body,  of  welded*  aluminium  alloy  conr 
strvction,  is  fitted  with  fixed  fins  at  the  rear  to  improve  the  aerodynamic 
characteristics  of  the  vehicle.  To  minimise  the  wei^t  penalty,  however, 
the  fins  are  not  large  enou^  to  make  the  airframe  aerodynamically  stable 
throxi^out  flight.  At  take-off  a  typical  vehicle  weighs  about  12,800  lb,  of 
which  propellents  ocoprise  approximately  11,400  lb.  The  Initial  acceleration 
is  thus  1*25g  (0,25g  excess)  and  the  velocity  and  altitude  at  all-burnt,  after 
145  sec  of  powered  fli^t,  are  about  11,500  ft  per  sec  and  400,000  ft  respec¬ 
tively,  During  the  ballistic  phase  peak  altitudes  of  about  500  miles  are 
attained. 

20 

Guidance  is  by  radio  command  link.  Deviations  from  the  desired 
f ligjit  path,  in  the  form  of  displacement  and  velocity  components  in  two 
orthogonal  planes,  corresponding  to  the  pitch  and  yaw  planes  wlien  the  vehicle 
is  ooirrectly  roll  stabilised,  are  derived  from  radar  or  optical  sources  and 
di^layed  on  cathode  ray  tubes  to  two  operators.  The  latter  are  able  to  vary 
the  vehicle  heading  to  keep  it  on  oovtrse  by  transmitting,  wiien  necessary, 
discrete  correction  signals  via  a  radio  link.  In  the  early  part  of  flij^t, 
when  the  desired  trajectory  is  vertical,  guidance  information  is  obtained 
from  two  M2  optical  trackers  situated  respectively  8000  yards  behind  and  to 
the  left  of  the  launcher.  Above  20,000  ft  the  desir*ed  flight  path  is  a 
strai^t  line  inclined  at  1°  50*  to  the  vertical  in  order  that  re*^nfery  may 
take  place  over  the  impact  area  50  to  60  miles  down  range.  For  this  part  of 
the  powered  ascent  phase,  guidance  information  is  obtained  either  from  a 
modified  S,G,R,584  radar,  with  its  aerial  dish  fixed  at  a  Q,E,  of  88°  10*,  or 
from  an  optical  telesccpe  system  in  which  the  rocket  motor  flame  is  tracked 
manually,  A  transponder  is  carried  in  one  of  the  vehicle  fin  pods  to  enable 
sufficient  radar  range  with  good  signal  to  noise  ratio  to  be  obtained. 

Filters  with  1,2  and  4  sec  time  constants,  which  can  be  inseirted  as  required 
by  the  operators  between  the  radar  resolvers  and  the  displays,  are  provided  to 
smooth  the  error  and  rate  signals. 

Cceanand  signals  are  transmitted  over  the  79  Ito/s  radio  link  in  the  form 
of  short  pulses  of  audio  frequency  modulation.  Four  a\adio  frequency  tones, 
which  correspond  to  left,  right,  vp  and  down  conmand  signals  respectively,  are 
available.  These  tones  may  also  be  transmitted  continuously  together  to 
command  self  descruction  in  the  event  of  the  vehicle  behaving  in  a  dangerous 
manner.  Discrimination  between  tones  is  aococpllshed  by  means  of  tuned  filters 
in  the  receiver.  Each  command  pulse  is  converted  by  the  command  units  shown  in 


♦Prototype  design  was  riveted, 
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Fig*  2  into  a  demand  for  0*25°  cshange  of  heading  in  the  pitch  or  yaw  plane. 

Fig*  3  shows  the  oircviit  diagram  of  the  ccnomand  vmitsj  these  effectively 
integrate  the  heading  changes  ccoinanded  in  the  two  planes  of  control,  their 
output  voltages  at  ai^  time  being  proportional  to  the  algebraic  sum  of  the 
heading  changes  demanded  since  take-off.  The  relay  shown  is  used  to  short 
circuit  the  integrator  capacitors  until  a  few  seconds  before  take-off.  To 
limit  the  rate  of  turn  demanded  by  the  guidance  system  to  one  degree  per 
second,  the  rate  at  which  pulses  can  be  transmitted,  via  any  one  audio 
frequency  channel,  is  fixed  at  four  per  second.  The  operators  guide  the 
vehicle  by  controlling  the  duration  of  the  pulse  trains  transmitted  from 
time  to  time. 

A  schematic  block  diagram  of  the  control  system  is  shown  in  Pig.  2. 

Roll  error  and  heading  angle  signals  are  obtained  from  two  free  g3a:osoopes 
fitted  with  A*C*  angular  displacement  pick-offs”*  ^  on  each  gimbal  axis. 

Signals  from  the  inner  gimbal  pick-offs  are  used  for  pitch  and  yaw  control 
and  from  both  outer  gimbal  pick-offs  for  roll  stabilisation.  Any  change  of 
roll  datum  after  the  gyroscopes  are  uncaged  is  thus  eq\ial  to  the  average 
drift  of  both  gyros  about  their  outer  gimbal  axes.  The  vehicle  is  roll 
stabilised  by  feeding  signals  \vhich  are  functions  of  roll  error  and  roll  rate, 
derived  from  a  phase  advance  network  and  associated  amplifier,  to  the  four 
motor  servos*  For  pitch  and  yaw  control  the  demodulated  gyro  signals  are 
passed  through  a  phase  advance  netv/ork  and  amplifier  and,  in  a  parallel 
branch,  added  to  the  oonnand  unit  output  signals  by  means  of  a  summing 
amplifier.  The  outputs  from  the  phase  advance  and  summing  amplifiers  are 
passed  throu^  the  gain  programming  vinit  and  combined  with  the  roll  control 
signals  in  the  summing  networks  of  the  motor  servo  amplifiers.  The  gain  pro¬ 
gramming  unit  consists  essentially  of  four  motor  driven  potentiometers  vhich 
vary  the  lateral  control  feedback  loop  gadLns  in  a  predetermined  manner  as  the 
fli^t  proceeds. 

The  theoretioal  design  and  analysis  of  the  control  system  should  take 
into  account  all  the  factors  which  influence  the  stability  of  the  system. 

The  following  are  the  main  ones  vrfiich  have  been  considered. 

(i)  The  vehicle  airframe  is  aerodynamiically  mstable  in  pitch  and 
yaw  during  the  greater  part  of  the  powered  flight  period.  In  other 
TTcrds  its  centre  of  pressure  is  then  situated  ahead  of  its  centre 

of  gravity, 

(ii)  The  vehicle  has  a  flexible  structvtre  and  hence  the  control 
gyroscopes  sense  local  deflections  of  the  structure  in  addition  to  the 
desired  motion  of  the  vehicle  as  a  whole, 

(iii)  The  mass,  manents  of  inertia  etc  of  the  vehicle  and  the  aero- 
dynamiic  farces  and  manents  acting  on  it  change  greatly  during  the 
powered  flight  period, 

(iv)  Liquid  propellents  with  free  surfaces  form  a  large  proportion 
of  the  total  mass  and  give  rise  to  ”sloshing"  effects. 

In  general  these  factors  increase  the  difficulty  of  designing  a 
satisfactory  control  system*  For  instance,  w4ien  the  vehicle  is  aerodynamic- 
ally  vtns table,  the  heading  control  feedback  loops  can  only  maintain  stability 
if  the  gain  exceeds  a  certain  minimum  value  wMch  depends  on  the  aerodynamiic 
mcment  derivative.  Otherwise  any  distuirbance,  e.g*  a  gust  or  wind  shear, 
would  ca\ase  the  vehicle  to  overturn  since  the  disturbing  mcment  would  inOTease 
more  rapidly  than  the  restoring  moment*  Hence  the  condition  of  maximum  aero¬ 
dynamic  instability  is  one  of  the  critical  design  cases*  On  the  other  hand, 
if  the  loop  gain  is  too  hi^  the  system  will  be  oscillatory,  as  with  any  feed¬ 
back  control  system*  In  the  case  of  a  flexible  missile  the  value  of  gain  which 
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makes  the  system  oscillate  depends  not  only  on  the  response  of  the  control 
system  components  but  also  on  the  body  bending  and  torsional  modes,  location 
of  the  gyroscopes  etc. 

The  purpose  of  the  analysis  vAiich  follows  is  to  determine  -vAiat  response 
the  control  system  components  should  have  and  what  valiies  of  gain  should  be 
specified  to  ens\ire  that  the  system  will  remain  stable  under  all  oonditions 
likely  to  occur  in  fli^t.  The  analysis  is  based  mainly  on  the  original 
estimates  of  aerodynamic  and  other  vehicle  parameters  but  the  effect  of  later 
modifications  is  referred  to  Yhere  appropriate.  Linear  servo-mechanism  theory 
and  frequency  response  methods  are  used  throvi^out.  Quasi-static  stability 
only  is  investigated;  in  other  words  the  vehicle  parameters  such  as  mass, 
speed  etc  are  assvimed  constant  for  each  instant  of  fli^t  considered  althou^ 
in  reality  they  vary  continuously  with  time.  This  assunption  that  stability 
depends  on  iixstantaneous  values  of  parameters  and  not  on  their  rates  of  change, 
is  not  thought  likely  to  lead  to  grossly  erroneous  conclusions,  as  the  rates 
of  change  are  generally  relatively  small,  but  it  should  be  borne  in  mind 
especially  Trtien  very  lav  frequencies  (belOTy  say  0,5  c/s)  are  considered. 

Motor  sei*vo  non-linearities  are  allowed  for,  to  some  extent,  by  using  exgperi- 
mentally  detennined  frequenoy  response  cvirves  for  four  different  anplitudes 
of  input  signal.  As  the  vehicle  is  roll-stabilised  and  symmetr^al,  apart 
from  the  fin  pods,  cross-coupling  effects  are  neglected  and  the  stability  of 
the  heading  control  system  in  one  plane  (yaw)  only  is  investigated.  The  pitch 
plane  control  system  is  assumed  to  be  identical  since  the  aerodynamic  estimates 
used  are  based  on  the  original  configuration  liiiiich  had  no  pods.  These  were 
introduced  as  a  modification  to  avoid  unacceptable  losses  in  the  transponder 
aerial  feeder  cable.  Since  the  pods  actually  make  the  airframe  scmneirtiat  more 
stable  in  yaw  than  in  pitch,  their  effect  is  to  give  the  yaw  control  system  a 
greater  margin  of  stability  than  is  iniicated  by  the  analysis.  The  gyroscopes 
are  assumed  to  be  located  in  the  electronics  bay  at  approximately  28  ft  from 
the  base  datum.  This  position  was  chosen,  after  preliminary  consideration  of 
the  stability  problem,  mainly  for  reasons  of  convenience  of  installation  and 
because  it  was  thougjit  that  the  vibration  environment  would  be  less  severe 
than  for  locations  nearer  the  rocket  motors.  Since  it  ■was  found  that  a  stable 
system  could  be  obtained  for  this  gyroscope  position  alternative  locations  ha-ve 
not  been  seriously  considered. 

The  first  part  of  the  Report  is  concerned  ■wi'th  the  stabili'ty  of  the  head¬ 
ing  con-trol  system  -vriiich  may  be  represented,  for  this  purpose,  by  the  sinplified 
block  diagram  shown  in  Fig,  4.  In  designing  a  closed  loop  oon'trol  system  it  is 
necessary  to  obtain  an  expression  for  "bhe  open  loop,  or  feedback,  "transfer 
function.  To  this  end  it  is  convenient  to  regard  the  system  as  consisting  of 
an  arbi"trary  number  of  interconnected  blocks  the  transfer  functions  of  ■vrtiich 
oan  be  detennined  fairly  readily.  The  "transfer  fvinction  of  a  "transmission  pa-tii 
through  a  system  cooponent  is  defined  as  "the  ratio  of  "the  Laplace  "transform  of 
the  output  signal  to  the  Laulaoe  "transform  of  the  input  signal.  The  Laplace 
transform  of  a  function  f (0  of  "the  real  variable  t  is  gi"ven  by 


OO 


f(t)  dt 


(1) 


e 


o 


■where  p  is  a  ccmplex  variable  cr  +  Jw.  The  "transfer  function  of  any  block  or 
component  "will  be  deno"ted  by  Y^(p),  or  simply  as  in  Fig,  4,  and  "tdie  corre¬ 
sponding  frequenoy  response  function  by  Y^(jw),  the  suffix  indicating  the 
block  concerned.  Transfer  functions  are  also  expressed  in  the  form 
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Yjj(p)  =  vdiere  is  a  gain  factor  iriiich  is  independent  of  frequency 

and  Gr^(p)  car  G-^Cjw)  deixstes  the  frequency  dependent  part.  Having  determined 

the  transfer  functions  of  the  individual  blocks,  the  open  loop  transfer  func¬ 
tion  is  obtained  sinply  by  adding  or  multiplying  these  as  indicated  by  the 
block  diagram.  This  iii5)lies  that  no  appreciable  interaction,  viiich  would 
alter  their  transfer  functions,  occurs  between  blocks  when  they  are  connected 
together  such  as,  for  exaii5)le,  loading  of  the  output  circuit  of  one  by  the 
ir5)ut  circtiit  of  the  next.  This  has  been  taken  into  account  in  drawing  the 
block  diagram  and  in  the  design  of  the  ccnponents.  Section  2  is  devoted  to 
the  missile  dynamics  transfer  function  YgCp),  "v^iich  relates  the  yaw  g3n’o- 

scope  pick-off  output  to  movements  of  the  yaw  motors  and  section  3  deals 
briefly  with  the  response  characteristics  of  the  rocket  motor  servos  and 
gyroscope  signal  demodulators.  In  sections  4,  5  and  6  the  stability  of  the 
heading  control  system  is  examined,  first  at  the  critical  time  of  maximum 
aerodynamic  instability  and  then  at  other  tiroes.  Propellent  sloshing  effects 
and  the  design  of  the  roll  stabilisation  ^stem  are  dealt  with  in  sections  7 
and  8  and  finally,  in  section  a  brief  account  is  given  of  ground  test  and 
flight  trial  resiilts, 

2  DERIVATION  OF  THE  MISSILE  DYNAMICS  TRANSBER  ITOCTION  IN  THE  YAW  PLANE 

2 

It  has  been  shoivn  by  Graham  and  Rodrignez  that,  for  the  purpose  of 
oalcxilating  the  forces  and  moments  due  to  small  oscillations  of  an  inconpres- 
sible,  invisoid  liquid  in  a  vertical  rectangular  tank,  an  eqiii valent 
mechanical  system  may  be  used.  This  consists  of  a  fixed  mass  plvis  an  infinite 
set  of  moving  masses  attached  to  the  tank  walls  by  means  of  springs  and  con¬ 
strained  to  move  horizontally.  Each  spring-mass  combinabion  represents  one  of 
the  characteristic  modes  of  oscillation  or  sloshing.  The  analysis  of  Ref*  2 
has  been  extended  to  the  case  of  a  vertical  cylindrical  tank  by  Staffoa:^  and 
numerical  val\ies  for  the  masses,  spring  rates  etc  for  Black  Kni^t  have  been 
calcxilated^.  In  the  case  of  Black  Kni^t  it  is  only  necessary  to  consider  the 
fundamental  sloshing  mode  for  each  tank  since  the  moving  masses  associated  with 
the  overtone  modes  are  small  enough  to  be  negligible.  For  instance  in  the 
early  part  of  fli^t  the  moving  masses  associated  with  the  ftjndamental  H.T.P, 
and  kerosene  sloshing  modes  are  13*0  and  7*4  slugs,  respectively,  whereas  the 
corresponding  second  mode  masses  are  only  0,4  0,2  slugs.  Fig,5  shows 

diagrasmatically  the  equivalent  missile  based  on  this  mechanical  analogue.  As 
the  propellents  are  consumed  during  fli^t  the  locations  and  magnitudes  of  the 
masses  etc  vary  with  time.  The  curves  of  Pig* 6  show  the  estimated  variation 
of  the  "effective”  mass,  lateral  moment  of  inertia  and  centre  of  gravity 
position  compared  with  corresponding  valxies  calculated  on  the  assxmption  that, 
at  any  instant,  the  propellents  remaining  in  the  tanks  are  solidified.  The 
"effective"  val\ies  are  calculated  fran  the  mass  etc  of  the  empty  vehicle  plus 
the  "fixed"  part  of  both  propellents.  Sprung  masses  are  not  included;  thxis 
the  difference  between  the  total  mass  and  the  "effective"  mass  at  any  time  is 
eqtial  to  total  of  the  sloshing  mode  moving  masses.  The  ounres  of  Pig*  6  are 
based  on  original  estimate b4  of  the  vehicle  parameters,  rates  of  propellent 
constmption  etc. 

The  complete  missile  dynamics  transfer  function,  including  aerodynamic 
body  flexure  and  propellent  sloshing  effects  all  together,  is  complicated  and 
inconvenient  to  use.  To  show  more  clearly  their  influence  on  the  stability 
of  the  control  system  these  effects  have  been  treated  separately  as  far  as 
possible.  It  has  been  assumed,  for  exanple,  that  body  flexure  effects  are 
independent  of  aerodynamic  and  sloshing  forces  and  depend  only  on  rocket 
motor  thrust  and  inertia  forces.  Similarly  propellent  sloshing  effects  are 
ignored  at  this  stage,  though  "effective"  values  of  centre  of  gravity 
position,  mass  etc,  as  defined  above,  have  been  \ised  thr cutout  in  calculating 
numerical  values  of  aerodyniunlo  parameters  and  transfer  function  coefficients. 
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Neglecting  propellent  sloshing,  the  equations  of  motion  in  yarr  of  a  rigid 
missile  may  be  Turitten  in  the  form 


m(v  +  xir)  = 


Y^v  +  T^r  +  ~  6  +  mg  t 


(2) 


Ct  =  N^v  +  N^r  -  H  6 


(5) 


where  the  dot  denotes  differentiation  with  respect  to  time, 

t  =  the  heading  angle  of  the  missile  with  respect  to  a  fixed 
vertical  datum; 

r  =  ijf  =  the  rate  of  turn  of  the  missile; 

6  =  the  deflection  of  the  two  rocket  motors  used  for  control  in 

■yaw  (see  Pigs.1  and  5)>  5  and  t  are  assumed  small  such  that 
sin  6  ^  b  and  cos  6^1  etc; 

T  ~  where  is  the  total  thrxist  of  the  motors, 

=  the  ccn^xinent  of  thrust  of  four  xmdeflected  rocket  motors  along 
the  missile  axis; 

u,  V  =  respectively  the  forward  and  lateral  ccnponents  of  velocity 
(see  Pig. 5); 


m,  C 


respectively  the  mass  and  lateral  moment  of  inertia  of  the 
missile; 

the  distance  of  the  centre  of  gravity  ft’om  the  motor  pivot 
station; 


g 


=  the  acceleration  due  to  gravity; 


Y^,  Y  =  the  aerodsmamic  normal  force  derivatives  with  respect  to 
^  V  and  r 

N^,  N^  =  the  corresponding  aerodynamic  moment  derivatives. 


Expressing  equations  (2)  and  (3)  in  terms  of  Laplaoe  transfonns,  assuming 
all  the  above  parameters,  except  5  and  v  and  their  derivatives^  to  be  constant, 
and  eliminating  v  yields  the  rigid  missile  d3mamics  transfer  function 


I 

6 


20 


lY  ) 


(4) 


where  i|f,  6  and  v  denote  Laplace  transforms  of  tjf,  S  and  v. 
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The  aerodynamic  derivatives  N  ,  N  and  Y  for  the  original  Black 

V  r  V  c 

Khi^t  configuration  at  zero  incidence  have  been  estimated  by  Stafford  and 
their  variation  Tdth  time,  for  the  assumed  standard  trajectory,  is  shown  in 
Pig*7»  Y^  is  asstaned  to  be  negligible.  It  is  evident  that  these  derivatives 

vary  considerably  during  the  powered  fli^t  period,  the  airframe  being 
vinstable  (N^  negative)  after  the  first  24  sec.  In  calculating  the  aero¬ 
dynamic  yawing  monents  about  the  C.G,  the  "effective”  centre  of  gravity 
(Pig,  6)  has  been  used  as  this  gives  a  greater  yawing  manent  than  would  have 
been  obtained  if  tiie  C.G.  based  on  instantaneously  solidified  propellents 
had  been  used.  The  resifLts  of  vdnd  tunnel  tests®  to  measttre  normal  force  and 
pitching  mcment  of  the  Black  Knight  configuration  at  various  valres  of  Mach 
number  and  Incidence  became  available  after  the  control  system  had  been 
designed  on  the  basis  of  Stafford’s  estimates.  As  the  tunnel  tests  indicated 
that  the  airframe  vras  less  unstable  than  had  been  estimated  it  was  decided 
that  no  alterations  to  the  system  were  necessary.  Pig. 8  shows  the  estimated 
variation  with  time  of  total  motor  thrust,  missile  acceleration,  velocity  and 
altitude  assxmed  when  the  aerodynamic  derivatives  of  Pig. 7  were  calculated. 
Numerical  values  for  the  coefficients  of  expression  (4)  can  be  conputed  as 
required  from  the  data  given  in  Pigs.  6,  7  and  8,  Per  exanple,  the  variation 
with  time  of  the  coefficient  T^/2Gf  which  gives  the  angular  acceleration  per 
unit  deflection  of  the  yaw  motOTs  in  the  absence  of  aerod3mamic  forces,  is 
shovm  in  Pig.  8  also. 

Since  the  missile  is  flexible  the  angle  detected  by  the  yaw  gyroscope 

Includes  a  conponent  tp  proportional  to  the  angular  deflection  of  the  body  at 

the  gyroscope  station  in  addition  to  the  desired  rigid  body  heading  angle 
Near  the  body  resonance  freqviencies  the  unwanted  feedback  due  to  this  body 
flexure  conponent  may  be  stifficlent  to  catise  control  system  instability  and 
precautions  may  be  necessary  to  prevent  this.  Thus  we  have 


=  ’I'  +  ’I'p  (5) 

where  it-  is  a  function  of  the  vehicle  mass  and  stiffness  distribution,  the 
£ 

location  of  the  g3rrosoopp,  the  lateral  force  P  and  reaction  couple  L  at  the 
motor  pivot.  It  is  assumed  that  is  independent  of  ilf  and  that  aerod3mamio 

forces  do  not  contribute  materially  tovTards  exciting  oscillations  at  the  body 
resonance  frequencies.  Por  small  5  the  exciting  force  and  coiple  are  given  by 


P  =  4t5+2M  d5 
^  c 

(6) 

and 

L  =  2  I  5 
c 

(7) 

where  I  and  M  =  respectively  the  manent  of  inertia  about  its  pivot  and  the 
°  °  mass  of  one  cocibustlon  chamber. 

d  =  the  distance  of  the  combustion  chamber  centre  of  mass  behind 
the  pivot  axis. 

In  the  case  of  Black  Kni^t  the  combustion  diambers  are  pivoted  about  axes 
udiich  pftss  very  close  to  their  centres  of  mass  (d  =!=  O)  and  hence  the  term 
2  d  6  in  equation  (6)  can  be  neglected. 
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It  is  assvimed  that  the  Laplace  transform  of  the  body  flexure  angle 
i|fp  is  given  by 


♦p 


H  „  P  +  H  -  L 
nP  nL 

— TT 

—  p  ♦  —  p 

n  0) 
n 


1  + 


(8) 


■where  w 
n 


2t^  ■vjhere  f  is  "the  nth  mode  resonance  frequency 
n  n  ^ 


n 


the  damping  factor  associated  with  -the  n-th  mode 


and  =  constants  for  "the  nth  mode  equal  to  2^^  x  the  angular 

deflection  airpli'tude  (in  radians)  at  "the  gyro  station 
at  resonance  per  lb  lateral  exciting  force  and  per  lb  ft 
couple,  respectively,  applied  at  "the  mobor  pivot  axis. 


Combining  equations  (6),  (7)  and  (8)  we  obtain  -the  body  flexure  part 
of  •the  missile  d3mamics  -transfer  function* 


^  V  "c 

TT 

+  72P 

n 


1  + 


(9) 


Bending  mode  shapes,  values  of  the  resonance  frequencies  f^  and  angular 

an5>lltudes  of  missile  deflection  at  the  gyro  station  due  to  lateral  forces  and 
couples  at  -these  frequencies  applied  at  the  mo-tor  pivots  ha-ve  been  calculated 
by  Donno®.  The  numerical  valxoes  are  given  in  Table  1  and  Jig*  9  shows  -the 
estimated  variation  -with  time  of  the  first  four  resonance  frequencies.  Par 
•these  calculations  the  missile  -was  assumed  -to  be  an  elastic  free-free  beam 
and,  in  the  absence  of  more  precise  information,  -the  danping  factor  was 

taken  as  0,01  (i.e.  of  critical)  for  all  modes.  Danping  of  -vibrating 
s-tructures  is  usually  assumed  to  be  hysteretic,  the  danping  force  being  pro¬ 
portional  to  displacement  and  in  j^ase  wi-th  velocity.  As  It  is  not  feasible 
to  ejpress  this  in  -transfer  function  form,  viscous  danping  has  been  assumed. 

It  has  been  shown  by  Bishop^  that,  viien  hysteretic  danping  is  small,  -the 
frequency  response  of  a  single  degree  of  freedom  system  is  -very  simileu?  to 
that  obtained  on  -the  assunption  of  viscous  danping. 

Equation  (9)  may  be  -written  in  -the  form 


n 
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where  =  ■?  T  and  the  frequency  of  motor 

^nL 

oscillation  at  which  excitation  of  the  nth  mode  by  motor  thrust  forces  is 
cancelled  by  that  due  to  the  inertia  reaction  coi5)le  L,  The  magnitude  and 
sign  of  h^  depend  on  the  mode  shape  and  the  position  of  the  gyroscope  rela¬ 
tive  to  the  nodes,  h^  is  positive  if  positive  motor  deflection  6,  shovm  in 
J^g«5*  tends  to  produce  positive  tj,  for  the  mode  considered.  Since  the 

proposed  gyro  position  is  forward  of  the  foremost  node  for  the  first  four 
modes,  is  positive  for  the  first  and  third  and  negative  for  the  second 

eind  fourth  modes.  Curves  are  given  in  Pigs.  10  and  11  which  show  how  the 
factors  h^  and  hg  and  the  frequencies  u^/27C  and  for  the  first  and 

second  bending  modes  vary  during  the  powered  fli^t  period.  The  higher 
modes  have  been  omitted  as  they  have  been  found  to  have  negligible  influence. 
These  curves  are  based  on  the  data  given  in  Table  1,  the  tharust/time  cvirve  of 
Pig. 8  and  an  estimated  combustion  chamber  moment  of  inertia  of  0.7*  slaig  ft^. 
General  frequency  re^onse  curves  corresponding  to  the  functions  (l  +  pvu^) 

•I 

and  •  y  — ■  ...I.-.,..  fQj.  ^  =0,01  are  given  in  Pigs,  12  and  13. 

.  Si  12  "" 

1  +— P  +"2^ 
n  or 
n 

Specific  transfer  functions  for  particular  times  and  the  corresponding 
frequency  re^onse  curves  can  be  obtained  as  required  from  equation  (lo)  by 
svibstituting  the  appropriate  values  of  h^,  u^  and  from  the  data  given  in 

the  figures.  Exanples  are  given  in  sections  4  and  5. 

By  equation  (5),  the  general  form  of  the  missile  dynamics  transfer 
function  Yg,  neglecting  propellent  sloshing  effects,  is  given  by  the  sum  of 

expressions  (4)  and  (10),  Thus 


3 


•  •  •  • 


(11) 


CHARAIC7EERISTICS  OP  THE  MOTOR  SERVOS  AND  GYROSCOPE  DEMQDDLATCRS 


Before  the  stability  of  the  heading  control  ^stem  as  a  vhole  is  considered 
a  brief  description  of  these  cotponents  is  necessary.  The  motor  servos  are 
designed  to  operate  at  a  hydravilic  pressure  of  3000  lb  per  sq  in,  and  to  have  a 
stall  torque  of  450  lb  ft.  This  latter  figure  was  chosen  to  make  liie  servo 
performance  insensitive  to  trimnion  friction  and  thrust  misalignment  torques, 
vAiich  at  the  time  were  not  acciirately  known,  and  to  meet  a  frequency  response 
specification  that  the  servos  should  be  capable  of  moving  the  motors  approxi¬ 
mately  sinusoidally  at  a  frequency  of  10  c/s  with  an  anplitude  of  Hydraulic 

power  is  obtained  in  fli^t  from  two  oil  punps,  driven  by  two  of  the  four 
prcpellent  punp  txirbines,  in  conjunction  with  an  hydraulic  accumulator.  When 


♦Actual  value  0,63  slug  ft 
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the  latter  is  fvilly  charged  the  pressure  rises  to  3»300p*s*i«  and  the  punrp 
delivery  is  by-passed  to  the  reservoir  through  an  off-loading  valve.  Oil  for 
the  servos  is  then  st^plied  from  the  accumulator  until  the  pressvtre  drops  to 
2,700  p»s*i#  •when  recharging  of  the  aocxmaulator  commences. 

The  layout  of  the  hydraulic  servo  is  shown  diagramnatically  in  Pig*  14* 

Oil  flow  -to  the  cylinders  is  controlled  by  a  single  stage  5/32"  diameter 
spool  valve  connected  by  a  torqiie  arm  to  the  shaift  of  an  E.12  torque  motor 
(commonly  known  as  "the  Laws  relay)”*  0  •which  is  totally  enclosed  and  imuersed  in 
oil.  The  ■val've  stem  profile  and  the  valve  liner  are  shaped  to  reduce  the  flow 
reaction  forces  to  a  low  level  and  -bhe  torque  motor  is  fitted  •with  a  restoring 
spring  to  gi've  optimum  loading  of  -this  oon^jonent  and  to  make  its  performance 
insensiti-ve  to  variations  in  pressure  and  flow  conditions'*"*  »”* 2,  The  actuator 
cylinders  are  connected  by  a  low  leakage  bleed  pipe,  not  shown  in  PLg,14»  to 
minimise  peak  press^ures  and  improve  servo  stability.  Motion  of  the  actuator 
pistons  is  •transmitted  to  the  rocket  motor  combus^tion  chamber  via  a  torque  arm 
and  two  s^truts  •with  ball  joints  at  each  end.  Backlash  and  •wear  are  automati¬ 
cally  taken  iip  wi-th  -this  arrangement. 

13 

AngtiLar  position  feedback  is  deri'ved  from  an  A.O.  pick-off  "the  rotor 
of  vAiich  is  •turned  through  •twice  the  combvistion  chamber  angle  by  a  mechanical 
l^mkage.  The  pick-off  ou'tput  is  demodulated  and  fed  back  to  one  of  the  seirvo 
anplifier  input  sunsning  resistors  as  indicated  in  Pig* 2.  Circuit  diagrams  of 
■the  demodulator  and  servo  amplifier  are  shown  in  Pigs*  15  end  l6.  These  com¬ 
ponents  and  "the  56  •volt  1,300  o/s  inverter,  •vshich  supplies  -the  pick-offs  and 
demodulators,  are  located  in  "the  electronics  bay  near  -the  front  of  -the  •vehicle. 
Primary  eleo-trical  power  supplies  are  obtained  from  batteries.  The  gain  of 
each  servo  demodulator  is  set,  in  conjunction  "with  its  associated  pick— off,  to 
give  1  volt  (d.G.  )  per  degree  of  motor  deflection  and  -the  gain  of  -the  servo 
aniplifiers  is  adjusted  to  give  a  differential  current  of  I6  mA  throuf^  "the 
torque  motor  coils  per  -volt  iigut.  Thus  a  D.C.  signal  to  any  of  "the  servo 
amplifier  1  MfJ  summing  resistors  results  in  a  steady  motor  deflection  of  one 
degree  per  volt  input.  Hence  =  1  degree  per  volt. 

Since  -the  dynamic  beha^viovir  of  •the  motor  servos  is  rather  nonrlinear,  a 
describing  function  based  on  experimentally  determined  frequency  response 
charaoteristios  of  a  prototype  unit  has  been  used  for  sys'tem  design  purposes 
instead  of  an  approximate  analytical  expression  for  •the  servo  -transfer  function 
Y^.  This  describing  function  G^(jw)  is  illus-trated  in  Pig*17  as  a  family  of 

gain  and  phase  curves  for  constant  anpli-tude  sinusoidal  ii^jut  signals  of  0w5» 
1,2  and  5  volts  wi'th  frequencies  ftom  1  to  100  c/s.  At  low  frequencies  -these 
iz^ts  give  motor  deflection  amplitudes  of  approximately  0.5,  1,  2  and  5  degrees 
respectively.  At  hi^er  frequencies  -vihen  -the  motor  angle  -waveform  is  not  sinu¬ 
soidal  -the  gain  and  phase  curves  are  based  on  -the  fundamental  cooponent  of  the 
output  waveform.  It  can  be  seen  from  Fig*  17  that  for  freqziencies  below  7  o/s 
•the  jhase  lag  increases  as  the  anpli-tude  decreases.  This  is  due  to  the  zise  of 
circular  ports  in  -the  valve  liner  and  to  valve  lap  and  stiotion  effects. 
Acceleration  and  velooity  sa-turation  are  responsible  for  -the  iixxreased  a-ttenua- 
tion  and  phase  lag  -vihich  are  apparent  for  -the  larger  amplitude  signals  at 
higher  frequencies.  Reduction  of  -the  feedback  gain  from  -thie  ipecified  figure 
of  16  mA  per  degree  of  motor  error  has  -the  effeot  of  incxreasing  -the  phase  lag, 
especially  for  small  anpli-tude  low  frequency  iiput  signals. 

The  gyroscope  demodulator  circuit  is  similar  to  -that  of  -the  servo  democu- 
lator  shown  in  Pig*15.  The  -transformer  ratios  are  different,  however,  since 
the  gain  is  hi^er.  The  gain  of  -the  pitch  and  yaw  demodulators  is  se-t,  in 
conjunction  -with  -their  respective  gyroscopes,  to  1*75  volts  per  degree  of  piok- 
off  rotation  and  -that  of  the  roll  demodula-tor,  -which  is  fed  fraa  -two  pick-offs 
in  parallel,  to  2*0  volts  per  degree.  To  a-void  inaccuracy  in  gain  se-tting  due 
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to  gyro  drift  effects  the  gyroscopes  are  oscillated  at  a  low  frequency 
with  known  angvilar  anplitude  v/hile  the  gain  potentiometer  RV2  (Pig#15;  is 
adjusted  vuatil  the  correct  demodulator  output  voltage  amplitude  is  obtained. 
The  demodulators  have  a  transfer  function  of  the  form 


'I'G  ^ 

Mhere  is  the  gain  in  volts  per  degree  and  is  a  time  constant  between 

1  and  2  milliseconds.  The  larger  value  has  been  used  in  this  analysis  and 

Pig*  18  shows  the  corre^onding  frequency  response  curves  for  G_(jw, 

5 

4  STABILITY  OF  THE  HEADING  CONTROL  SYSTEM  AT  THE  TIME  OF  MAXIMUM 

AERaDYNAtaC  INSTABILITY 

The  stability  of  the  heading  control  system  may  now  be  examined  using 
the  transfer  and  describing  functions  from  the  previous  sections.  The  object 
is  to  establish  the  required  gains  and  frequency  characteristics  of  the 
remaining  control  system  components  shown  in  the  block  diagram  of  Pig#4»  the 
phase  advance  and  simming  amplifiers.  For  a  system  of  this  type,  the  most 
convenient  method  to  ^ase  for  this  purpose  is  the  Nyquist  diagram  and  stability 
criterion"*^.  The  latter  enables  the  number  of  roots  with  positive  real  parts 
in  the  chareioteristio  eqmtion  of  a  closed  loop  control  system  to  be  deter¬ 
mined  from  the  open  loop  frequency  response  and  transfer  functions.  The 
I^yquist  diagram  is  a  complex-plane  plot,  or  vector  locus,  of  the  open  loop 
transfer  fijnction,  as  the  conplex  variable  p  is  allowed  to  assume  values  on  a 
oontovcr  which  encloses  the  entire  right-hand  half  of  the  p-plane.  If  N  is  the 
number  of  counter-clockwise  encirclements  of  the  point  (-1,j0)  by  this  locus 
the  I^yquist  stability  criterion  states  that 


N  =  Z  ~  P  (13) 


•viiiere  P  =  the  nmber  of  poles  of  the  open  loop  transfer  function  in  the 
right  half  of  the  p-plane  (l*e,  with  positive  real  parts) 

and  Z  =  the  number  of  roots  of  the  closed  loop  characteristic  equation 
with  positive  real  parts. 

For  a  stable  system  Z  must  be  zero.  Inspection  of  the  Nyquist  diagram  also 
gives  useful  information  about  the  relative  stability,  or  damping,  of  the 
closed  loop  system  and  about  vAiat  Increase  or  decrease  of  loop  gain,  for 
example,  would  be  permissible.  This  margin  of  stability  can  be  expressed 
numerically  in  terms  of  the  gain  and  phase  margins  associated  with  specific 
valiaes  of  loop  gain  etc.  If  these  margins  are  reasonably  large  (say  6  db 
and  30^)  the  closed  loop  resp>onse  of  the  system  will  normally  be  adeqviately 
damped. 


Before  a  Nyquist  diagram  can  be  drawn  it  is  necessary  to  derive  an 

expression  for  the  open  loop  transfer  function.  As  the  required  phase 

advance  and  svniiming  amplifier  transfer  functions  Y,  and  Y^  are  unknown  at 

4  5 


this  stage,  the  assiiiption  is  made  that  they  have  the  simple  forms  I^p 
,  respectively,  vhere  and  are  gain  factors  to  be  determined. 


and 

This 


is  equivalent  to  assTmning  that  the  deflection  of  the  rocket  motors  is  pro- 
pxartional  to  heading  error  plus  rate  of  change  of  heading  gyro  oulput.  To 
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avoid  making  the  first  Nyquist  diagram  too  ccmplicated,  the  hody  flexure  terms 
in  the  missile  dynamics  transfer  function  (11)  have  been  emitted.  The  rigid 
missile  transfer  function  at  t  *  80  sec,  iiriien  uN  ,  tlie  aerodynamio  moment  per 
vinit  incidence,  has  its  maximum  value,  is  given  oy 


8,36p  0.696 


p^  +  0,3l6p^  -  9»45p 


0,16 


(14) 


This  expression,  v/hich  is  derived  from  eqtiation  (4)  and  the  data  given  in 
Figs«6  to  8,  has  two  poles  in  the  ri^t  hand  half  of  the  p^lane.  The 
cotrre^onding  open  loop  transfer  function  for  t  =  80  sec  is 


^lV3(\ 


I^y^Cp)  (I^P  K^)  (8.36p  0,696) 

(1  +  0.002p)  (p^  +  0.3l6p^  «  9,45p  +  0.16) 


(15) 


where  G^(p)  is  defined  by  the  family  of  gain  and  phase  cvtrves  shown  in  Pig*17« 

Curve  A  of  Pig,  19  sho\7s  the  Nyquist  diagram  corresponding  to  expression  (15) 
for  the  1°  amplitude  motor  servo  frequency  response  G^  ( j<o)  and  the  arbitrary 

gain  values  of  KjK^K^  =  1  and  =  4«  These  assumed  gain  factors  are 

equivalent  to  4  degrees  of  motor  deflection  per  degree  of  heading  error  and 
1  degree  of  motor  deflection  per  degree  per  second  rate  detected  by  the  gyro. 
From  the  ccuplete  Nyqxiist  diagram,  shown  as  inset  A  in  Pig,  19,  it  can  be  seen 
that  the  vector  locus  encircles  the  point  (-1,o0)  twice  in  the  cloclcvise 
direction.  Thtis  N  =  -2  and.  since  P  =  2  for  the  open  loop  transfer  fxnKxbion 
(15),  Z  =  0  by  equation  (135*  Hence  for  the  rigid  body  case  the  system  is 
stable  with  these  values  of  gain.  It  can  be  seen  from  curve  A  that  the  phase 
margin  is  47°  at  1.3  o/s  and  the  gain  margins  are  -10,5  db  at  0^107  c/s  and 
22,4  db  at  17«5  o/s.  This  implies  that  if  the  error  loop  geiin  were  reduced 
by  more  than  10,5  db,  i,e,  a  factor  of  3«35,  the  system  would  become  unstable 
due  to  aerodynamic  moment  increasing  more  r^idly  than  the  restoring  moment 
from  thrust  deflection.  Alternatively  if  the  gain  were  increased  by  more  than 
22,4  db,  a  feiotor  of  13*2,  the  system  would  become  unstable  in  an  oscillatory 
manner. 

Ourve  B  of  Pig,  19  and  the  corresponding  inset  show  the  effect  of  neglect¬ 
ing  the  gravity  term  gN^/C  in  equation  (4),  The  open  loop  transfer  function 

then  becomes,  for  the  same  gain  factors 


G-j(p)(P  +  4)  (8.36p  +  0,696) 
p(1  +  0.002p)(p^  +  0.31 6p  -  9*45) 


(16) 


This  function  has  one  pole  in  the  ri^t  half  of  the  ccnplex  plane  (P  =  I) 
and  one  at  the  origin.  The  conplete  Nyquist  diagram,  inset  B,  shows  one  nega¬ 
tive  encirclement  of  (-1,j0)  and  hence  Z  =  0  which  indicates  that  the  system  is 
stable,  as  would  be  expected.  The  gain  and  phase  margins  are  practically  the 
same  as  for  loctas  A  since  the  loci  differ  appreciably  only  for  frequenoies 
lower  than  0*1  c/s.  Locus  B  tends  to  00  /~270^  as  w  tends  to  zero  from  the 
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positive  direction.  In  plotting  the  loci  A  and  B  a  logarithmic  (dh)  scale 
has  been  used  for  magnitudes  greater  than  2  and  a  linear  scale  for  smaller 
magnitudes.  This  has  been  done  to  enable  the  Iovt  frequency  part  to  be  shown 
more  fully  without  using  an  unduly  small  linear  scale.  It  does  not  affect 
the  number  of  encirclements  of  nor  the  gain  and  phase  margins,  thou^ 

it  does  distort  the  shape  of  the  loci. 

The  effect  of  including  the  first  two  body  flexure  modes  is  illustrated 
by  Pig, 20  which  shows  the  Nyquist  diagram  for  t  =  80  sec  corresponding  to  the 
open  loop  transfer  function 


Y^YgY^CY^+Y^) 


(p)  (p  +  4) 

(1  +  0,(X)^) 


a,36p+  0.696 


p^  +  0.3l6p^  -  9.45P  +  0.16 


6,15  X 10"^ 

p" 

1  + . . p 

(2vtx28,0)  . 

{2%  X  10.5) 


0.02p 

+  t;  Vo  v  h-1 


2  "  2%  x  10.5 


1,25  X  10“^ 

1  +  9 

L  (27C  X  24.0)  . 

'] 


0,02p 

+  V"  V  +  1 


(27tx29i3)^  2x  X  29.3 


.... (17) 


This  expression  is  derived  from  eqxjations  (II  )  and  (15)  and  the  data 
given  in  Pigs.  6  to  11,  The  same  arbitrary  values  of  gain,  4  and  1  seo“'',  have 
been  used  and  the  1°  amplitude  case  of  Pig. 17  has  been  used  for  G^Cjw).  In 

this  case  the  system  is  unstable  since  the  open  loop  transfer  function  has 
two  poles  in  the  ri^t  half  plane  as  before  but  the  net  number  of  encircle¬ 
ments  of  the  point  (-1,j0)  is  aero.  The  number  of  aeros  is  thertfore  2,  by 
equation  (13)|  indicating  that  the  closed  loop  characteristic  equation  would 
have  two  roots  7rf.th  positive  real  parts. 

One  way  of  stabilising  the  system  is  to  change  the  simple  functions 
initially  assumed  for  Y^  and  Y^  to  obtain  considerably  more  phase  lag  at  the 

fundamental  body  flexure  frequency.  Examination  of  the  effects  on  system 
stability  of  a  number  of  possible  filters  indicated  that  the  insertion  of  a 
third  order  low  pass  filter  with  a  break  frequency  of  3  o/s  vrauld  give 
acceptable  stability  throughout  flight.  Prequency  response  curves  of  such 
a  filter  are  shown  in  Pig,  21 .  The  open  loop  transfer  function  with  this 
filter  included  becomes 


^1  V3^\  +  Y5) 


(p  +  4)  G^(p) 


(1  +  0.002p)  (1  +  O.OlSpy’ 


[Missile  dynamics  transfer  function  as  in  (17)  ] 

....(18) 
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Cvirve  A  of  Pig.  22  shov/s  part  of  the  Njrquist  diagram,  oorr emending  to  this 
function,  plotted  for  frequencies  above  0*01  c/s.  The  lower  frequency  part 
is  the  same  as  that  of  Fig.  20.  It  is  evident  from  this  diagram  that  the 
system  is  now  stable  since  the  couple te  Nyquist  diagram  encircles  the  point 
(-■Ijjo)  twice  in  the  negative  sense,  as  for  the  rigid  body  oase  shown  in 
ilg.19  (curve  A),  Cotparison  with  the  latter,  however,  indicates  that 
insertion  of  the  fd.lter  has  resulted  in  substantial  reduction  of  the  phase 
and  gain  margins.  The  phase  margin  is  now  29°  at  1.25  c/s  and  the  gain 
margins  are  7.1  db  at  3.25  o/s  and  -10,7  db  at  0,12  o/s 

Since  at  body  resonance  frequencies  the  gain  of  the  phase  advance 
acplifier  is  much  larger  than  that  of  the  sunming  amplifier,  it  is  not 
necessary  to  include  the  filter  in  the  latter  and,  in  practice,  it  is  simpler 
and  more  convenient  to  insert  it  in  the  phase  advance  amplifier  only.  The 
open  loop  transfer  function  then  becomes 


0.  (p)  r  -|  ^ssile  dynamics  ■) 

....(19) 

Locus  B  of  Pig,  22  is  the  Nyqtiist  diagram  for  this  oase,  ©le  system  is 
stable  with  a  phase  margin  of  32°  at  1,45  c/s  and  gain  margins  of  6,0  db  at 
3.5  c/s  and  -10.7  db  at  0,12  c/s.  The  transfer  function  of  the  actual 

phase  advance  amplifier  differs  slightly  from  that  incorporated  in  expression 
(19)  since  A.C.  coupling  is  used  throu^out  to  avoid  saturation  by  the  D.C. 
oonponent  of  the  gyro  demodulator  signal  and  to  minimise  amplifier  drift 
effects.  The  actual  open  loop  transfer  function  is 


^1^3^^  +  ^5) 


(p) 

(1  +  0,00^) 


4  + 


16,8 


—  (1  +  2^,1p)^  (1  +  0,0l8p)^. 


ell 


r  Missile  dynamios 
transfer  function  [ 
as  in  (17) 


.... (20) 


The  2f.1  sec  time  constants  associated  with  the  interstage  coupling  networks 
shown  in  the  circuit  of  Pig,  23  affect  the  gain  and  phase  at  low  frequencies 
as  Indicated  in  Pig,  21f.  Curves  A  correspond  to  the  transfer  fuiwtlon  in 

expression  (19)  and  curves  B  to  that  in  (20),  Locus  C  of  Pig. 22  shows  how 
the  Nyquist  diagram  is  affected.  The  phase  margin  is  new  34°  at  1,45  o/s  and 
the  gain  margin  at  3.5  c/s  is  unchanged  at  6,0  db  but  the  low  frequency  gain 
margin  is  1,1  db  less  at  -9.6  db  at  0,11  o/s.  The  latter  would  be  still 
further  reduced  if  the  coupling  network  time  constants  were  made  smaller.  As 
they  stand  those  values  of  gain  and  phase  margin  are  large  enou^  to  assure 
satisfactory  closed  loop  stability  and  adequate  damping  but  further  redixction 
would  be  undesirable, 

Conparing  the  Ifyquist  diagrams  of  Pigs, 20  and  22,  it  is  evident  that 
introduction  of  the  stabilising  filter  has  restilted  in  a  considerable  reduction 
in  the  magnitude  of  feedback  at  the  second  body  flexure  mode  fSrequenoy  of 
29,3  o/s.  The  gain  margin  associated  with  .this  mode  is  about  I6  db  whi^  should 
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be  more  than  adequate  even  if  the  strvictural  danping  factor  'vrere  lovTer  than 
the  assvuned  value  of  1^.  A  lovrer  fundamental  mode  dan^jing  factor  Trould  not 
affect  the  stability  of  the  system  since  the  phase  angle  at  10»5  o/s  is 
such  that  feedback  is  negative  at  this  frequency.  The  influence  of  hi^er 
bending  modes  on  the  stability  of  the  system  appears  to  be  negligible. 

Finally  the  effect  of  motor  servo  non-linearities  must  be  considered. 

The  Nyquist  diagrams  shovm  in  Pigs.  19,  20  and  22  are  based  on  servo  fre¬ 
quency  response  measurements  for  1°  aiiplitude.  Par  demanded  ainplitude, 
hofrever,  the  measured  phase  lag  at  1.5  and  3  o/s  is  4°  and  5°  greater  than 
for  the  1°  case  though  the  gain  is  imchanged.  This  extra  phase  lag  redijoes 
the  phase  margin  to  30^  at  1.45  c/s  and  the  gain  margin  to  5«2  db  at  3*2  c/s. 
It  is  clear  from  Pig.  17  that  at  low  frequencies  the  servo  phase  lag  increases 
as  the  input  signal  anplitude  decreases.  If  the  system  stability  margins  were 
inadequate  this  could  result  in  small  amplitude  limit  cycle  hunting.  Par 
irput  signal  amplitudes  larger  than  1  volt  (i.e.  1°  demand)  the  1cm  frequency 
gain  and  phase  margins  77111  be  greater  than  those  indicated  in  Pig.  22  but  at 
hi^er  frequencies  the  attenuation  and  jiiase  lag  are  greater.  For  fjistance, 
at  the  fundamental  flexure  frequency  of  10.5  o/s  the  phase  lag  is  5^°  greater 
for  5°  than  for  1°  demanded  ampLitude.  It  can  be  seen  from  the  Nyquist 
diagram  that  this  extra  phase  lag  does  not  make  the  proposed  system  unstable 
nor  the  stability  margin  inadequate.  In  any  case  it  is  very  unlikely  that 
such  large  amplitudes  and  rates  of  motor  movement  v/iU  ever  be  demanded  in 
normal  flight.  It  is  therefore  concluded  that  the  proposed  system  will  be 
satisfactory  for  the  conditions  expected  to  occur  at  the  time  of  maximum 
aerodynamic  instability. 

5  HEADING  CCNTRQL  SYSTEM  STABILITY  AT  OTHER  TIMES 


Practical  considerations  make  it  desirable  that  control  system  changes 
in  fli^t  should  be  restricted  to  variations  of  the  gain  factors  ^  and  K^; 

changes  of  filter  time  constants  etc  would  involve  greater  complexity  and 
hence  potentially  decreased  reliability*  It  will  now  be  shown  that  the 
system  arrived  at  in  the  previous  section  is  satisfactory  from  this  point 
of  view. 

Stability  immediately  before  all-burnt  (t  =  145  sec)  and  immediately 
after  clearing  the  launcher  (t  =  o)  are  considered  first  as  these  are  the 
times  at  which  the  vehicle  mass,  inertia  etc  have  their  minimum  and  maximum 
values.  Since  the  dsmamic  pressure  is  very  small  for  both  these  cases,  aero¬ 
dynamic  forces  and  moments  have  negligible  influence  on  control  system 
stability.  Thus  the  rigid  body  transfer  function  (4)  reduces  to 


I  =  -  (21) 

8  2Gp^ 


••2 

where  T^2C  =  16.7  sec  at  t  =  145  and  5.85  at  t  =  0  as  shown  in  Pig.8,  Since 
the  angular  acceleration  of  the  rigid  missile  per  unit  motor  deflection  at 
t  =  145  sec  is  approximately  double  that  at  t  =  80  sec  it  is  necessary  to 
reduce  the  internal  control  system  gain  by  a  factor  of  two  to  obtain  the  same 
margin  of  stability.  As  the  latter  is  not  excessive  at  t  =  80  sec  it  was 
decided  to  increase  the  gain  and  phase  margin  at  all-bvcrnt  by  reducing  K  and 


K5  ty 


a  factor  of  four  since. 


in  the  absence  of  aerodynamic  moments  there  is 


no  sharply  defined  lower  limit  to  the  gain  idiich  can  be  used.  Too  lew  a 
value  of  gain  would,  hosTever,  make  the  response  of  the  vehicle  to  guidance 
command  signals  sluggish  and  poorly  danped  and  this  would  adversely  affect 
guidance  accuracy.  Pig.  25  shows  the  Nyquist  diagram  for  t  =  145  sec,  plotted 
from  the  open  loop  transfer  function 
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where  (p)  denotes  liie  motor  sei^  transfer  function  for  1  °  aii5)litude. 

Higher  body  flexure  modes  are  omitted  as  they  do  not  affect  stability 
significantly.  It  is  evident  from  Fig. 25  that  the  system  is  stable  with  a 
phase  margin  of  38°  at  Cu9  c/s  and  a  gain  marg^Ln  of  14.8  db  at  3*4  o/s.  If 
the  internal  gain  were  kept  the  same  as  at  t  =  80  sec  the  phase  margin  would 
be  10^  and  the  gain  margin  2.8  db  which  would  be  quite  inadequate.  There 
appears  to  be  no  danger  of  instability  at  the  fundamental  flexure  frequency 
even  if  the  strijctural  dan^jing  should  be  less  than  has  been  assumed. 

The  %qTaist  diagram  for  conditions  immediately  after  take-off,  shown 
in  Hg.  26,  has  been  plotted  for  the  same  values  of  internal  gain  as  those 
selected  for  t  =  80  sec  ft:*om  the  open  loop  transfer  function 


lV3(\  +  V  =  “ 


G/p) 


4  + 


I6.8p" 


1  +  0.002^  [_  ^  4,1p)^(l  +  0.0l8p)^ 
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It  can  be  seen  that  the  system  is  stable  v/ith  a  phase  margin  of  37°  at  1*1  o/s 
and  gain  margins  of  7*8  and  7*5  ^  at  3*3  and  I6.5  c/s  respectively.  Feedback 
at  the  fundamental  body  flexure  frequency  is  predominantly  negative  and  at  the 
first  overtone  frequency  the  phase  angle  is  such  that  neither  motor  servo  nonr* 
linearity  nor  a  moderate  reduction  in  structural  damping  will  make  the  system 
unstable.  Feedback  at  the  higher  mode  frequencies  is  negligible. 

From  the  three  oases  considered  so  far  it  appears  likely  that  stable 
fli^t  TiTould  be  obtained  if  the  internal  gain  factors  and  were 

kept  constant  at  4  and  1  respectively  until  shortly  after  t  =  80  sec  and  were 
then  gradually  redvjced  to  1  and  0.25  at  145  sec.  In  practice,  however,  lower 
values  of  gain  at  take-off  are  preferable  and  it  can  be  seen  from  Fig. 26  that 
the  gain  at  t  =  0  could  be  reduced  by  5  db  without  any  redviction  in  phase 
margin.  Such  a  redTiction  would  have  the  effect  of  reducing  the  motor  deflec¬ 
tions  at  launch  due  to  gyro  and  anplifier  drifts  and,  in  the  event  of  the  gain 
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prograimie  motor  failing  to  rvm  after  take-off,  the  lower  gain  wotild  be  a 
better  conpromise  for  a  fixed  gain  system  and  ■would  give  stable  fli^t 
throughout.  Immedia'tely  after  release  a  transient  disturbance,  invol'ving 
lateral  motion  of  the  vehicle  rear,  is  likely  to  occur  due  to  'the  cccibined 
effects  of  •wind,  thrust  ineqvialities  and  misalignments,  motor  deflections 
at  take-off  etc.  This  is  unimportant  provrided  that  -the  lateral  movement  is 
less  "than  the  available  laimcher  clearance  and,  to  ensure  that  'this  require¬ 
ment  is  met,  limit  svrf.'tches  are  incorporated  viiioh  pre'vent  'the  ■vehicle  being 
released  if  ■the  deflection  of  any  motor  exceeds  Reducing  ■the  gain  at 
take-off  ■thus  reduces  the  probabili'ty  of  this  limit  being  reached.  The  gain 
programme  shown  in  Fig«27  Tvas  ■therefore  selected.  The  heading  error  gain 

increases  linearly  frcm  2«25  at  t  =  0  to  4  at  t  =  75,  remains  constant 

\mtil  t  =  85  and  then  decreases  vinifarmly  to  1  at  t  =  145  sec.  The  rate  gain 
varies  similarly  being  al'ways  one  quarter  of  ■the  above  figures.  Due 


to  the  characteristics  of  -the  linkage  connecting  ■the  potentiometer  ■wipers  to 
the  dri'ving  shaft  •the  actual  gain  programme  differs  slightly  from  "that  shown 
in  Pig,  27. 


Al'thou^  it  is  fairly  obvious  from  ■the  ■vmy  in  ■which  ■the  aerodynamic  and 
O'ther  parameters  vary  ■with  time  -that,  for  a  rigid  vehicle,  ■the  proposed  gain 
programne  v/ill  give  a  s'table  and  well  damped  system  throughout  ■the  powered 
fli^t  period  it  is  not  immediately  ob'vicus  that  instability  due  to  body 
flexure  will  be  a-voided.  To  shovr  this,  vector  loci  have  been  plotted  in 
Pig, 28  ■«?hich  gi-ve  the  variation  ■wi'th  time  of  ■the  miagnitude  and  phase  of  ■the 
feedback  at  body  flexure  frequencies.  Loci  and  Ag  are  dravm  for  the  funda¬ 
mental  mode  for  Jhe  1°  and  5°  servo  ampli-tvide  cases  respectively  and  B  for  the 
first  overtone  1  servo  case.  As  before  a  s'tructural  damping  factor  of  of 
critical  has  been  ass'umed.  To  find  out  if  'the  system  is  xmstable  due  to  'the 
body  flexure  mode  considered  at  any  particular  timae  it  is  necessary  to 
determine  the  approximate  shape  of  ■the  Nyquist  diagram  for  frequencies  close 
to  the  body  resonance.  This  can  be  done  by  simply  drawing  a  circle  ■with  the 
■vector  for  the  time  of  in'terest  as  dlameter31.  Par  example,  circle  gives 

the  approximate  shape  of  the  liPyquist  diagram  near  frequency  f^  at  t  =  1 09  sec 
for  the  5°  ampli'tude  case.  The  system  is  stable  provided  tha-t  this  circle 
does  not  enclose  the  point  (-Ijdo).  This  procedure  does  not  give  -the  exact 
shape  of  ■the  Nyquist  diagram  but,  as  can  be  seen  from  ■the  second  mode  in 
Pig,26  far  example,  it  gives  reasonable  accuracy  near  ■the  resonance  frequency 
■sdien  ■the  damping  is  law  and  when  feedback  due  to  body  flexure  is  large  cam^‘ 
pared  ■with  that  d\ie  -to  rigid  body  motion.  It  can  be  seen  from  curve  A-j  that 
for  ■the  1°  ampli'tude  case  •the  system  is  stable  and  would  remain  so  even  if 
■the  s'tructural  damping  factor  ■were  considerably  lower  than  'the  assxmned  value. 
Par  the  5°  case  locus  Ag  shows  ■that  ■the  system  is  stable  but  that  if  the 

damping  factor  ■were  less  ■than  half  the  assumed  ■value  instability  mi^t  occur 
be'tween  100  and  115  sec.  Howe'ver,  since  sudi  large  ampll'tudes  or  rates  of 
motor  movement  and  such  small  ■values  of  s-tructural  damping  are  ■very  unlikely, 
the  risk  of  instabili'ty  is  -very  smiall.  The  ■worst  case  for  ■the  second  mode 
(curve  B)  appears  to  be  t  =  0  -.»hich  has  already  been  ^ecifically  dealt  ■with 
(Pig,26)  and  shown  to  be  stable.  It  is  therefore  concluded  ■that  ■with  the 
system  proposed  instability  due  to  body  flexure  should  not  occur. 

6  THE  EPPECT  OP  CHANGES  IN  THE  BODT  BENDING  MODES 


After  ■the  basic  theoretical  design  of  the  lateral  con'trol  system  had 
been  completed,  as  outlined  in  -the  previous  sections,  sUbs'tantlal  alterations 
to  the  ■vehicle.  s'tructTxre  ■lirere  in'trodxjced.  These  invol'ved,  inter  alia, 
replacement  of  ■the  original  riveted  body  ■wi'th  stringers  by  a  hea'vier  gaxage 
■welded  ■version  ■without  s'tringers.  It  ■was  therefore  necessary  to  re-assess 
the  effects  of  body  flexure  on  the  control  system  on  the  basis  of  re'visod 
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estimates  of  bending  mode  frequencies  etc.  Table  2  gives  the  revised  figures, 
vrtiich  are  reprodvioed  ftrom  Table  3  of  Ref.  15.  Ccnparison  iTith  Table  1  shows 
that  in  the  early  part  of  flight  the  bending  ft*equeix3ies  have  increased  but 
towards  all-burnt  they  are  considerably  lower.  The  farced  vibration  anrplitudes 
at  the  gyro  station  are  generally  somewhat  less  than  far  the  original  structure. 

Pig* 29  shows  the  Nyquist  diagram  for  the  tanks  full  case  (t  =  O)  based 
on  these  revised  estimates,  the  other  elements  of  the  control  system  being  the 
same  as  before.  As  it  was  thotight  that  the  VTelded  vehicle  might  have  lower  • 
strvictural  danping  than  the  riveted  one  a  damping  footer  of  0*3?o  ^ras  assumed 
for  this  case.  It  is  clear  from  the  Nyquist  diagram  that  even  with  this  low 
value  of  damping  the  system  is  stable  with  an  anple  margin.  In  particular  it 
may  be  observed  that  the  Increase  in  the  first  oveiixjne  frequency  has  resulted 
in  a  substantially  increased  gain  margin  for  this  mode  (c.f.  Pig.  26). 

To  show  that  the  system  is  stable  throughout  fli^t,  loci  similar  to 
those  of  Pig,  28  have  been  plotted  in  Pig, 30  xising  the  revised  estimates  for 
the  fundamental  bending  mode,  A  strvotural  danping  factor  of  0,9^  has  been 
assumed.  These  loci,  the  interpretation  of  vsbich  has  been  disciissed  in  the 
previous  section.  Indicate  that  the  system  should  remain  stable  throughout 
the  powered  fli^t  period.  Plence  it  was  concluded  that  no  alteration  to  the 
proposed  control  system  was  necjessary  on  account  of  the  improved  structure, 

7  THE  EPPECTS  OP  ERCEELLENT  SLOSHING 


Expressions  for  the  transfer  function  relating  rigid  body  heading  angle 
to  rocket  motor  deflection  with  sloshing  effects  iixsluded  have  been  derived 
for  two-tank  ballistic  missiles  by  Smith"*®  and  Ward"* 7  on  the  basis  of  the 
mass-spring  mechanical  analogue  referred  to  in  section  2.  To  keep  the  order 
of  the  transfer  function  low  only  the  fundamental  sloshing  mode  for  each  tank 
has  been  incltided;  the  effect  of  the  overtone  modes  is  considered  negligible. 
In  the  absence  of  aerodynamic  farces  the  transfer  function  has  the  form 
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=  re^ectively  the  mass  and  spring  rate  associated  -with  the 
fundamental  sloshing  mode  for  the  H.  T*P.  tank 

kg  =  the  mass  and  spring  rate  for  the  fundamental  kerosene  tank 
mode 

■^2  -  distances  of  masses  m.  and  m^  behind  the  effective  centre 

of  gravity 

R^  =  mass  and  radius  of  gyration  of  eroty  missile  plus  "fixed" 
part  of  both  propellents.  (M  R^  =  effective  lateral 
moment  of  inertia  C)  °  ° 

a  =  oonponent  of  acceleration  along  missile  axis  (assumed  vertical). 

The  variation  with  time  of  the  fundamental  H.T.P.  and  kerosene  sloshing 
frequencies  fp  eind  f^^  is  shown  in  Plg,31,  fran  Ref ,4,  "vAiere 


In  the  early  part  of  fll^t  fp  and  fj^  are  nearly  equal,  varying  from 

1,13  c/s  at  t  =  0  to  1«3  c/s  at  60  sec.  Thereafter  they  diverge  slowly, 
fp  rising  to  a  peak  of  2»8  o/s  at  135  sec  whereas  fj^  reaches  its  maximum  of 

1*9  o/s  at  120  sec.  Towards  all-burnt  fp  and  fj^  fall  rapidly  towards  zero. 
The  transfer  function  of  expression  (2^)  can  be  factorised  to  give 
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•srfiere  Kg  =  1^/20  as  before  and  oo^,  c*^,  etc  are  respectively  27t  x  frequencies 
f^,  f^  and  f^,  f^  at  which  this  function  has  zeros  and  poles.  The  variation 
of  these  with  time  from  launch  is  shown  in  Pig,  32.  It  can  be  seen  that  f 


and  f  are  nearly  equal  throughout  powered  fli^t  and  vary  in  a  similar 
c 

manner  to  fp  vdille  f^  and  f^  follow  closely  the  variation  of  fj^.  Since  in 

deriving  the  spring-mass  analogue  the  propellents  have  been  assumed 
inviscid,  the  transfer  function  (26)  becomes  infinite  at  frequencies  f  and 
f^  and  the  Nyquist  diagram  consists  of  tvTo  infinite  semi-circles.  It  San 
be  seen  from  tiie  shape  of  the  Nyquist  diagrams  obtained'^?  that,  in  the 
absence  of  propellent  daaping,  the  stability  of  the  system  depends  on  the 
order  of  frequencies  f  ,  f,  ,  f  and  f ,.  Referrii^  to  Pig, 32  there  are 
three  cases s-  ^ 


(1) 

(2) 

(5) 


f^<f^<f^<f^  (i,e,  zero,  pole,  zero, 
t  =  0  to  85  sec  and  113  to  140  sec 

fl^<f^<f^<f^  (i.e.  zero,  pole,  pole, 
t  =  85  to  113  sec 

f^  <  f^  <  f^  <  f^  (i,e,  pole,  zero,  zero, 
t  «  140  to  145  sec. 


pole). 

System  stable. 

zero). 

System  unstable 

pede). 

System  unstable 
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To  obtain  quantitative  information  about  the  damping  of  the  sloshing 
modes  under  closed  loop  conditions  it  is  necessary  to  evaluate  the  roots  of 
the  characteristic  equation  derived  from  the  closed  loop  transfer  function. 
Since  the  fuiviamental  sloshing  frequencies  are  rather  lovr,  it  was  thought 
that  the  lags  associated  with  the  motor  servo,  filter,  gyro  demodiiLator  etc 
ootild  be  neglected  in  the  first  instance.  Then,  to  confirm  this  and  to 
obtain  seme  indication  of  the  effects,  if  any,  of  these  lags,  a  simple 
second  order  transfer  function,  having  about  the  same  total  gain  and  phase 
over  the  sloshing  frequency  band,  was  incltBied  and  the  roots  again  evaluated. 

Neglecting  conponent  lags  the  open  loop  transfer  function  is 


Y 

o 


(27) 


and  the  characteristic  equation,  from  the  dencminator  of  the  closed  loop 
transfer  function  1^1  -  Y^) ,  is 


P^(P^  +  <^)  (P^  + 


0  . 


(28) 


The  roots  of  eqviation  (28)  have  been  obtained,  by  means  of  a  digital  OGcputer, 
for  parameter  valites  at  intervals  of  ten  seconds  or  less,  assuming  the  para¬ 
meters  to  be  constant  for  each  instant  considered.  The  gain  factors  Kg, 

K^KjK^  and  K^K_K_  were  assumed  to  vary  as  shown  in  Pigs.  8  and  27  and  the 


frequencies  f^,  f^,  f^  and  f^  as  in  Pig. 32.  Six  complex  roots,  corresponding 

to  three  oscillatory  modes  with  frequencies  f  ,  f  and  f  and  damping  facto(rs 

X  y  z 

^  and  were  obtained  for  each  instant  considered.  Pig»33  shows  how 


these  frequencies  and  damping  factors  vary  as  the  propellents  are  used  vp. 

^x*  i^tural  frequency  of  the  heading  control  system  in  the  absence 

of  aerodynamics  and  conponent  lags,  increases  from  0*57  o/s  at  take-off  to 
0*9  c/s  at  80  sec  diminishing  again  to  0*56  o/s  at  145  sec  due  to  the  programmed 
gain  reduction.  The  corresponding  damping  factors  are  0*44,  0*6?  and  0*51» 

The  other  two  are  closed  loop  sloshing  modes  and  it  can  be  seen  that  f  and  f 

Jr  ^ 

vary  in  a  similar  manner  to  fj^  and  fp  respectively.  The  damping  factor  of 


what  may  be  loosely  called  the  closed  loop  kerosene  sloshing  mode  starts  very 
small  (0,1^  of  critical)  and  positive,  increases  to  0*027  at  100  sec  and  then 
falls  to  zeiro  at  140  sec.  During  the  last  five  seconds  it  is  negative,  as  v/ould 
be  expected  since  f^  ^  very  small  (f  0,003)  v/hich  implies  a  slowly 


diverging  oscillation.  The  H.T.F.  closed  loop  sloshing  mode  damping  factor  ^ 
is  positive  and  greater  than  0,01  vmtil  66  sec.  As  would  be  expected  there  is  a 
period  of  instability  from  83  to  113  sec  but  it  can  be  seen  from  Pig#33  that 
the  greatest  negative  damping  factor  is  0,005  at  100  sec.  Prom  113  sec  to  all- 
burnt  this  mod.e  is  stable. 
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The  effect  of  including  control  system  cooponent  lags  is  illustrated 
in  Pig* 34*  In  this  case  the  open  loop  transfer  ftinction  Yra.s  taken  as 


K^KgK^d^p  -f  K3)  |-(p^  4  og)  (p^  +  (cg) 


P^(1  +  0.04B7P  +  0,0019^^)  +  <^)  (p^  + 


•] 


(29) 


■SThich  gives  an  eighth  order  characteristic  eqmtion.  The  additional  second 
order  term  in  expression  (29)  was  chosen  to  represent  the  effective  lag  of 
the  control  system  ccanponents  as  it  has  about  the  same  total  gain  and  phase 
lag  for  frequencies  fran  1  to  3  c/s.  CcjEparison  of  Pigs. 33  and  34  shovTs 
that  ■with  this  extra  lag  included  the  closed  loop  sloshing  frequencies  fy 


and  f^  are  practically  the  same  as  f^  and  f^  in  the  early  part  of  flight 
but  are  slightly  higher  after  about  1CX)  sec.  The  damping  factor  of  the 
kerosene  mode  varies  in  a  similar  manner  to  S  and  is  generally  somewhat 

y 

larger.  The  period  of  instability  associated  v/ith  the  H.T.P.  mode  occurs 
sli^tly  later,  between  90  and  117  sec  and  the  maximum  negative  value  of 
^  (-0*006)  occurs  at  110  sec.  Between  120  and  140  sec  ^  is  appreciably 

less  than  and  a  second  very  short  period  of  instability  occurs  around 

135  sec  =  -0,0006)  \ihen  f^  has  its  maximum  valroe  of  3*27  o/s.  As  this 


is  the  main  difference  between  the  resvilts  obtained  frcm  the  sixth  and 
eighth  order  equations  it  is  reasonable  to  conclude  that  only  \Then  the 
sloshing  mode  frequency  exceeds  2,5  c/s  is  its  closed  loop  danping  factor 
affected  adversely  to  an  appreciable  extent  by  control  system  oonponent 
lags. 


TThen  the  closed  loop  sloshing  mode  damping  factor  remains  negative 
for  a  i)erlod  any  disturbance  v/ill  res\ilt  in  oscillations  of  increasing 
amplittide,  Por  instance  at  100  sec,  triien  =  -0*005,  oscillations  at 

2,1  o/s  would  build  up  to  double  their  initial  anplitude  in  10,5  sec,  if 
the  parameters  remained  constant.  The  angjlitude  would  continue  to  increase 
until  either  the  damping  factor  became  positive  or  scxne  non-linear  effect, 
such  as  wave -breaking,  caused  anplitude  limiting  to  occur.  Even  when  the 
sloshing  mode  damping  is  positive  it  is  evident  that  if  propellent  oscil¬ 
lations  were  excited,  by  aerodynamic  disturbances  or  guidance  conmands, 
these  would  take  a  long  time  to  decay  as  the  damping  factors  are  so  small 
especially  ^  in  the  early  part  of  fli^t.  In  view  of  this  it  was  decided 
that  the  guidance  command  repetition  frequency  should  be  changed  from  two 
per  second,  as  originally  proposed,  to  fo\ir,  since  periodic  disturbances 
at  or  near  one  of  the  sloshing  frequencies  might  cavise . forced  oscillations 
to  build  vq)  to  large  amplitude.  The  frequency  of  four  per  second  was 
chosen  as  it  is  above  the  highest  fundamental  sloshing  frequency  and  below 
the  lov/est  body  flexure  frequency. 

The  effects  of  aerodynamic  forces  and  of  propellent  viscosity  have 
so  far  been  ignored  in  the  discussion  of  sloshing.  Vard‘'7  has  shown  that, 
although  the  influence  of  aerodynamic  forces  on  the  frequency  and  damping 
of  the  closed  loop  sloshing  modes  is  small,  they  tend  to  make  the  system 
more  unstable  when  one  of  the  sloshing  modes  is  already  luistable  in  the 
absence  of  aerodynamics.  Sane  idea  of  the  effect  of  propellent  viscosity 
may  be  obtained  by  asstaning  each  spring-mass  to  be  associated  with  a 
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viscous  damper,  although  the  mechanical  analogy  is  strictly  applicable  only 
for  inviscid  fluid.  It  has  been  sho\-,'n  by  Smith that  if  at  100  sec  either 
the  H,T,P,  or  kerosene  fundamental  sloshing  mode  is  assumed  to  have  a  damp¬ 
ing  factor  of  0,05  associated  v/ith  it,  the  other  being  undan^jed,  the 
unstable  (z)  mode  becomes  positively  damped  \vith  a  daii5)ing  factor  of  +0#010 
or  +0*006.  Thus  it  appears  that,  yliile  viscous  damping  in  either  tanlc  is 
beneficial,  the  increase  in  closed  loop  danping  factor  is  less  than  that  of 
the  open  loop  mode  and  that  damping  in  the  H,T,P,  tank  is  slightly  more 
effective  for  this  particular  unstable  mode.  The  Black  Knight  tanks  are  not 
smooth-walled  but  contain  stiffening  rings  spaced  about  18"  apart  in  the 
H.T.P.  tanlc  and  12"  apart  in  the  kerosene  tanlc.  Experimental  measurements 
by  Lineham’^9  have  sho\7n  that  in  -vvater  svich  a  ring  gives  a  danping  factor  of 
0*02  vdien  at  the  surface  and  0*003  vdien  9"  belorr  it.  For  a  smooth  3  ft 
diameter  tank  the  danping  factor  v/as  found  to  be  0, 001 .  In  flight  the  danp- 
ing  v/oxild  be  less  than  this  due  to  the  effect  of  acceleration  but  this  would 
be  partly  counterbalanced  by  the  kinematic  viscosity  of  kerosene  being 
greater  than  that  of  v/ater. 

The  case  for  additional  anti-sloshing  devices  in  Black  Khi^t  is  thus 
rather  inconclusive  as  the  period  of  instability  is  fairly  short  and  the  rate 
of  divergence  of  sloshing  oscillations  appears  likely  to  be  slow.  However, 
to  ensure  the  greatest  chance  of  success  in  the  early  flight  trials,  it  vra.s 
decided  that  seven  equally  spaced  rows  of  fourteen  baffle  plates,  each  3” 
square,  should  be  welded  to  the  H.T.P.  tank  Trails  betvTeen  42"  and  79”  from 
the  lower  datum  to  increase  the  sloshing  mode  damping  between  88  and  II7  sec* 
Consideration  was  also  given  to  tiie  possibility  of  increasing  the  sloshing 
mode  damping  under  closed  loop  conditions  by  modifying  the  heading  control 
system  transfer  function  electrically  but  this  did  not  appear  to  be  practic¬ 
able  ^.'ithout  the  use  of  extra  feedback  transducers  such  as  lateral  accelero¬ 
meters.  Such  extra  conplicatlon  vras  considered  to  be  un;}ustifled.  After 
all,  even  if  heading  an^e  oscillations  of  moderate  amplitude  did  occxir  for 
a  period  during  the  ascent  they  w'ould  be  unlikely  to  result  in  catastrophic 
failure.  At  a  later  stage  in  the  trials  prograirme  it  may  be  considered 
desirable  to  flight  test  a  Black  loiight  without  anti-slosh  baffles,  which 
Involve  a  small  wralght  penalty,  to  determine  experimentally  if  they  can  be 
safely  left  out, 

8  THE  ROLL  STABILIS.\TICK  SISTEM 


If  large  deviations  from  the  correct  roll  attitude  were  to  occur  in 
flight,  due  to  aerodynamic  aix3  other  disturbing  torques  acting  on  the  vehicle, 
cross-cotroling  betvTeen  the  azimuth  and  elevation  guidance  channels  would  make 
accxirate  grxLdance  difficult.  It  has  been  shown^*^  that  a  total  phasing  error 
of  15°  can  be  tolerated  without  serious  effects  on  the  performance  of  the 
guidance  system  though  considerably  less  than  this  is  obviously  desirable. 
Allowing  a  somewhat  pessimistic  total  of  10°  for  gyro  drift  and  phasing  error 
in  the  guidance  radar  set,  it  appears  that  tp  to  5°  roll  error  with  respect 
to  the  gyro  datum  v/ould  be  acceptable.  To  minimise  the  effects  of  gyro  drift, 
signals  from  the  outer  gimbal  pick-offs  of  both  gyroscopes  are  used  to  feed  a 
common  demodulator.  In  addition  to  averaging  the  drift  of  the  two  gyros, 
this  arrangement  has  the  advantage  that  the  system  v/ould  still  function,  with 
reduced  gain,  if  the  roll  pick-off  of  one  gyro  became  inoperative  due  to  an 
open-circuit  faxilt. 

The  maximum  transient  disturbing  torque  due  to  motor  thrust  inequalities, 
fin  misalignments  etc  vras  originally  estimated  to  be  1000  lb  ft,  the  bulk  of 
this  being  due  to  aerodynamic  forces  at  the  time  of  maximum  dynamic  pressure. 
Since  the  motor  thrust  acts  at  a  radius  of  0*81  ft  the  control  moment  available 
in  roll  varies  fVom  226  lb  ft  per  degree  of  motor  deflection  at  sea  level  to 
256  above  100,000  ft.  Thus  the  maximum  motor  deflection  likely  to  be  required 
to  overcome  the  estimated  disturbing  torque  is  about  4°  and  hence,  to  keep  the 
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error  •within  5°>  a  control  system  stiffness  of  not  less  than  0*8°  of  motor 
deflection  per  degree  of  roll  error  is  desirable.  The  possibility  of  using 
integral  of  error  con'tar'ol  irras  considered  but  it  riaa  decided  "that  the 
possible  benefit  ■would  be  insufficient  to  justify  the  extra  oonplioatlon 
involved.  There  should  be  no  steady  state  rolling  mctoent  and  it  rr&s  felt 
that  an  integrator  -would  be  ano-fcher  potential  source  of  drift. 

In  analysing  -the  stability  of  the  roll  stabilisation  system  -the  proce¬ 
dure  adopted  is  similar  to  -that  used  for  "the  heading  con-trol  system.  The 
aim  is  to  determine  -what  -transfer  functions  should  be  inserted  be-tween  -the 
gyro  demodulator  emd  -the  rocket  motor  servos  to  ensure  -that  -the  system 
remains  stable  throughout  -the  povTered  fli^t  period.  The  problem  is  sinpler 
-than  for  -the  heading  con-trol  system  since  it  has  been  shown^^  that  -the  bulk 
of  the  propellents  do  not  rotate  -vrfien  the  vehicle  rolls.  It  may  therefore  be 
assumed  that  the  mcment  of  inertia  in  roll  remains  constant  at  -the  value  for 
■the  empty  -rehicle.  This  -was  estimated  to  be  57»2  slug  ft^.  The  maximum 
value  of  -the  aerodynamic  roll  dan^jing  coefficient  Lp  was  estimated  to  be 

75  ft  per  radian  per  sec,  at  50  sea  after  take-off,  v?hlle  -the  propellents 
oon-tribute  a  small  amount  of  -viscous  dan^jlng.  For  practical  purposes  this 
is  almost  negligible  since  -vdth  tanks  fvill  it  amounts  -to  about  3*5  It*  ft  per 
radian  per  sec. 

Considering  first  the  vehicle  as  a  rigid  body,  -the  equation  of  rolling 
motion  is 


A?  +  Lp  ^ 


(30) 


where  ^  =  roll  angle  wi-th  respect  to  roll  da-fcum 

A  =  moment  of  inertia  in  roll 
=  roll  daiiplng  coefficient 
R  =  radius  at  -iihich  motor  thrust  acta 

=  ccnponent  of  motor  deflection  tending  to  roll  -the  vehicle 
=  roll  disturbing  torque. 

The  rigid  missile  -transfer  function  relating  roll  angle  to  motor  deflection 
is  therefore 


£ 


TR 

A 


p  +XP 


(31) 


where  TIV'A  =  227,  22+2  and  257  sec  and  I^A  =  0*06,  1*35  sec**'  and  zero 
for  t  =  0,  50  and  >120  sec  respectively, 

■where  Kg  and 

rate  gain  factors  at  our  disposal  and  and  Y^  are  the  motor  servo  and  gyro 
demod\jlatar  -transfer  functions,  the  open  loop  -transfer  function  becomes 


are  roll  error  am 


Assuming  that  6^^  =  "^1^^%^  + 
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Y 

o 


Y,  K,  f-  (Kg  .  K^) 
(1+0,  CI02p)  +  A  P 


(32) 


assiming  the  gyro  (iemociulator  transfer  function  to  be  the  same  as  before. 

Pig* 35  shov/s  part  of  the  Nyqxoist  (iiagram,  plotteci  frcm  expression  (32),  for 
four  cases  -which  illus-trate  the  effects  of  the  variation  of  motor  servo 
response  v/ith  irput  anplitucie  ard  of  the  rigici  bo(3y  -transfer  function  with 
time.  In  all  four  cases  the  assumeci  gain  factors  are  K,K,K^  =  1*7  (degrees 

of  motor  deflection  per  degree  roll  error)  and  K^Iv^Kg  =  0,1  sec  -which  give 

the  maxlmvim  possible  stiffness  consistent  v/ith  an  adequate  margin  of  stability. 
Curves  A,  B  and  C  are  (iravm  for  -the  -g-®,  1°  and  5°  servo  anplitude  cases  -wi-th 
the  maximvan  value  of  TR/A  (i»e,  t  >  120  sec)  vrfiile  curve  D  shov/s  the  effect  of 
Lp  at  50  sec  for  the  amplitude  case,  Conparison  of  ctLTves  A  and  D  shows 

that  Lp  only  affects  the  shape  of  -the  Nyquist  diagram  significantly  at  fire- 

quenoies  below  2  c/s.  If  plotted  to  a  smaller  scale,  it  vrould  be  seen  that 
curve  D  becomes  asynptotic  to  the  negative  imaginary  axis  as  the  flrequency 
tends  to  zero.  In  all  four  cases  the  system  is  stable  -VTlth  adequate  gain  and 
phase  margins.  For  the  1°  anplitude  case  -the  gain  margin  is  10,1  db  and  the 
phase  margin  32°.  The  latter  is  reduced  to  25°  (ciirve  A)  or  27°  (curve  D)  for 
the  ampli-tude  case  but  the  gain  margin  is  increased  slightly.  For  the  5° 
anplitude  case  -the  phase  margin  is  43°  and  the  gain  margin  7*9  db. 

Since  torques  applied  to  the  vehicle  by  the  rocket  motors  -\/ill  produce 
torsional  deflections  of  the  structure,  it  is  necessary  to  investigate  to  vrhat 
extent  these  affect  the  stability  of  the  roll  stabilisation  system.  For  small 
motor  deflections  -the  applied  torque  is  TR6p  due  to  -tliriist  and  -the  conponent 

of  -torsional  coiple  due  to  angular  acceleration  of  the  combustion  chambers  is 
4Iq6jj  sin  5f°*  The  transfer  function  relating  roll  gyro  output  to  motor 

deflection  6^,  may  be  written  in  the  form 


2 

P 


+ 


J^(TR  +  4Iq  p^  sin  5|°) 


2^: 


2 


P 


(33) 


•where  the  first  term  is  the  rigid  body  -transfer  function  and  the  subsequent 
ones  the  effects  of  torsional  modes; 

=  2'Xf^^  -vThere  f^^  is  the  nth  mode  torsional  frequency 

=  the  danping  factor  associated  ivi-th  the  nth  torsional  mode 

=  a  constant  for  -the  nth  torsional  mode  equal  to  2^^  x  the 

angular  anplitude  of  torsional  deflection  (in  radians)  at 
resonance  at  -the  gyro  station  per  lb  ft  amplitude  of  -torsional 
covple  applied  at  the  motor  pivot  station. 
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The  fundamental  torsional  frequency  f_.  v/as  originally  estimated  to 

-6  ' 

be  50  c/s  vdth  =  15*6  x  10  assuming  to  be  ^fa  and  the  system  vras 

designed  on  this  basis.  Curve  A  of  Fig,  36  shovTs  the  Nyquist  diagram  includ¬ 
ing  this  mode;  it  is  directly  conparable  with  curve  A  of  Fig,35»  being  drawn 
for  the  same  assumed  values  of  gain  and  the  servo  an^slifier  case,  and 
corresponds  to  the  open  loop  transfer  function 


Y 

o 


(1.7  +  0.1p) 

R 

(1  +  0. 00^) 


4.6  X  10"^  (1  +  0.19  X  10"^ 

1 

p 


.  0,02 
1  +  vT;-  o'  P  + 


314.2 


(314.2)‘ 


•  •  •  • 


(34) 


The  feedback  amplitude  at  higher  torsional  mode  frequencies  was  assumed 
to  be  negligible.  It  is  evident  frran  the  Nyquist  diagram  that  the  system  is 
stable  as  the  point  (-1,j0)  is  not  encircled  and  the  margin  of  stability  near 
the  fundamental  torsional  frequency  is  considerable.  This  is  increased  for 
larger  servo  amplitudes  due  to  the  greater  servo  phase  lag  at  50  c/s. 
Referring  to  expression  (34)  >  the  term  (l  +  0, 19  x  10“^  p2)  is  zero  at 
36,6  c/s;  in  other  v/ords  at  this  frequency  the  torsional  excitation  due  to 
motor  thrust  is  cancelled  by  that  due  to  combustion  chamber  inertia  effects. 
At  lower  frequencies  thrust  forces  provide  the  main  excitation  and  at  higher 
frequencies,  such  as  we  are  concerned  Td.th  here,  inertial  effects  are  the 
most  inportant. 

To  avoid  troubles  due  to  noise,  especially  demodtilator  ripple,  it  is 
necessary  to  limit  the  bandwidth  of  the  phase  advance  amplifier  by  means  of 
low  pass  filters.  The  larger  the  time  constants  of  these  filters  the  smaller 
the  bandwidth  will  be  and  the  less  noisy  the  output  but  the  greater  the  extra 
phase  lag,  A  compromise  is  necessary  therefore  and  time  constants  of  3 
1.5  milliseconds  vrere  chosen  to  give  adequate  rejection  of  demodulator  ripple 
and  a  still  acceptable  phase  margin.  Curve  B  of  Pig,  36  sho\7s  the  effect  of 
these  filters  on  the  Nyquist  diagram  and  Pig,37  gives  the  frequency  response 
and  transfer  function  of  the  proposed  roll  phase  advance  amplifier.  The 
phase  margin  is  only  20°  for  the  amplitude  case,  thou^  for  larger  servo 
amplitudes  it  is  greater,  and  the  gain  margin  is  6,3  db.  If  larger  time 
constants  vrere  used  greater  noise  rejection  would  be  obtained  but  it  would 
be  necessary  to  reduce  the  stiffness  of  the  system.  It  was  concluded,  on 
the  basis  of  the  original  torsional  mode  estimates  that  a  maximum  stiffness 
of  1,7°  of  motor  deflection  per  degree  of  roll  error  cotild  be  used  if  the 
frequency  response  of  the  phase  advance  amplifier  \7as  as  shown  in  Pig, 37* 

23 

The  torsional  modes  were  subsequently  re-calculated  to  allow  for  the 
effects  of  fin  bending  and  of  changing  to  the  welded  teink  structure  referred 
to  in  section  6,  The  fundamental  torsional  mode  for  this  was  estimated  to  be 
63,4  c/s,  with  =  4*3  X  10"°,  and  the  fin  modes  to  be  39  and  41  c/s.  The 

effect  of  these  changes  on  the  stability  of  the  system  vreis  re-assessed  by 
calc\ilatlng  the  magnitude  and  phase  of  the  feedback  vector  at  the  resonance 
frequencies.  The  fin  modes  were  found  to  have  negligible  effect  and  the 
magnitude  of  the  feedback  vector  at  63.4  o/s  was  found  to  be  0,208  with  a 
phase  angle  of  -356°.  It  was  tlierefore  concluded  that  no  instability  3ue  to 
body  torsion  was  to  be  expected. 
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9  GROUND  TEST  AND  FLIGHT  TRIAL  RESULTS 

It  has  been  shown  in  the  preceding  sections  how  the  main  features  of 
the  control  system  were  arrived  at  by  using  linear  servo-mechanism  theory 
and  the  best  available  estimates  of  vehicle  parameters.  To  coii55lete  this 
account  it  is  relevant  to  Include  a  brief  simanary  of  the  control  system  tests 
subsequently  carried  out  and  to  indicate  their  influence  on  the  evolution  of 
the  system. 

7?hile  design  eind  manufacture  of  prototype  test  vehicles  and  conponents 
were  proceeding,  a  series  of  single-plane  fixed-peirameter  simulator  studies 
of  the  heading  control  and  roll-stabilisation  systems  was  carried  out  using 
an  electronic  analogue  conputer  and  a  prototype  hydraulic  servo.  In  these 
studies,  which  have  been  reported  in  detail  elsevhere^^ the  dynamic 
performance  of  both  systems  was  examined  for  a  number  of  points  on  the 
standard  trajectory  over  a  range  of  gain  factors  and  their  responses  to  a 
variety  Of  stimuli  were  recorded.  No  major  design  changes  were  found  to  be 
necessary  as  a  result  of  these  simulator  studies. 

At  a  later  stage  closed  loop  tests  of  the  control  system  were  carried 
out  using  conplete  vehicles  on  a  static  firing  test  stand  with  and  without 
the  motors  thrusting.  The  object  of  these  tests  was  to  determine  whether  any 
undesirable  effects,  due  to  structxiral  constraints  etc,  would  occur  if  the 
gyro  feedback  loops  v/ere  closed  prior  to  take-off.  In  the  first  series  of 
tests,  with  the  motors  not  thrusting,  the  heading  control  systems  were  found 
to  be  satisfactory  but,  when  the  roll-stabilisation  loop  vras  closed,  the 
motors  oscillated  with  small  amplitude  at  about  65  c/s.  In  an  attempt  to 
ascertain  the  cause  of  these,  open  loop  measurements  vrere  made  in  which  the 
rocket  motors  vrere  oscillated  and  the  resultant  motion  of  the  vehicle  in  the 
vicinity  of  the  gyros  was  measured  by  accelerometer  type  pick-ups.  Conpari- 
son  of  the  angular  anplitudes  indicated  by  the  latter  and  by  gyroscopes 
shovred  that  at  65  c/s  the  gyro  outputs  were  very  much  larger.  The  nutation 
frequency  of  the  gyros  then  in  use  (Reid  and  Sigrist  type  GS75  iik.l)  ^ras 
measured  and  foxind  to  be  about  65  c/s  with  a  damping  factor  of  less  than  0«5^» 
It  was  therefore  suggested  that  an  alternative  gyroscope  (Smiths  Series  7 
Mod. 4)  with  a  nutation  frequency  of  112  c/s  should  be  tried  in  the  hope  that 
this  wo\ild  cure  the  trouble.  In  the  first  aeries  of  tests  with  the  new  gyro¬ 
scope  both  heading  control  and  roll-stabilisation  systems  functioned  satis¬ 
factorily  in  the  absence  of  motor  thrust.  Subsequently  the  roll-stabilisation 
system  again  became  noisy,  the  rocket  motors  tending  to  buzz  at  about  40  o/s» 
illthough  this  vms  probably  caused  by  deterioration  of  the  gyro  performance 
under  vibration,  it  was  apparent  that  a  reduction  in  the  gain  of  the  roll  rate 
amplifier  at  frequenoies  between  20  and  1(X)  c/s  (see  Pig, 37)  v/ould  be  bene¬ 
ficial  provided  that  this  could  be  obtained  without  excessive  low  frequency 
phase  lag.  To  meet  these  requirements  a  tuned  L.C.  filter  giving  more  than 
20  db  attenuation  between  40  and  50  c/s  and  a  small  low  frequency  phase  lag 
v/as  designed  and  inserted  between  the  roll  phase  advance  and  motor  servo 
amplifiers.  The  small  increase  in  lor/  frequency  phase  lag  due  to  the  filter 
made  it  necessary  to  reduce  the  stiffness  of  the  roll-stabilisation  system 
(K^ICjKg)  to  0,8  to  maintain  adequate  gain  and  phase  margins.  As  mentioned 

in  section  8,  this  lower  stiffness  is  still  sufficient  to  give  acceptable 
performance  even  if  the  maximum  estimated  roll  disturbing  torque  occurs. 

Insertion  of  the  filter  was  found  to  give  satisfactory  control  system 
performance  in  the  test  stand  both  in  the  absence  of  thrust  and  dviring  static 
firing  of  the  motors.  During  pre-laurwh  testing  of  the  first  flight  vehicle 
(BKOI)  at  'Jocanera,  however,  audible  60  c/s  noise  was  found  to  occur  and  it 
was  decided  that  increased  attenuation  at  this  frequency  would  be  desirable. 

A  two-stage  L.C.  filter,  the  circuit  and  frequency  response  of  which  are 
shown  in  Pig,38  was  made  and  inserted.  At  first  the  trouble  seemed  to  have 
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"been  ciored  but  next  day  noise  again  occurred.  This  vas  found  to  be  due  to 
one  of  the  gyroscopes  having  beccme  noisy  and  after  this  had  been  replaced 
no  further  noise  trouble  was  experienced.  BK. 01  was  fired  on  7th  September, 
1958  and  the  control  i^stem  performed  satisfactorily  in  flight^^. 


Pig* 3 9  shoii/s  the  I’yquist  diagram  for  the  roll  stabilisation  system 
with  the  two-stage  filter  inclvided  and  the  stiffness  reduced  to  0,8.  It 
can  be  seen  that  the  system  is  stable  with  gain  and  phase  margins  of  5»0  db 
and  27°  respectively  for  the  -5°  servo  airplitude  case.  For  the  1°  case  these 
increase  to  o, 0  db  and  34°*  No  change  was  made  in  the  rate  gain  factor 

thovi^  this  could  well  have  been  reduced  by  3*5  db  without  appreciably 

affecting  the  dynamic  performance  of  the  system. 


To  check  the  bending  and  torsional  mode  estimates  a  series  of  full- 
scale  resonance  tests  has  been  carried  out  on  a  Black  ICnight  vehicle  with 
welded  tanks  (BIC, 02)30,  Tfith  tanks  fiill  the  meastcred  bending  mode  resonance 
frequencies^^^  were  generally  in  quite  good  agreement  with  the  estimates  but 
the  aaplitudes  at  the  gyro  station  vrere  lower.  The  fundamental  mode  (9«83  c/s) 
danping  factor  v/as  about  but  for  the  second  and  third  modes  values  of  1*^ 
and  3*1^  were  obtained.  These  damping  factors  \7ere  calculated  from  energy 
balance  considerations;  the  corresponding  figures  derived  from  polar  plots 
of  anplitude  and  phase  from  selected  pick-vps  were  0.9^,  1.1%  and  1#3^» 

V/ith  tanlcs  enpty  the  ftandamental  resonance  frequency  (23.5  c/s)  agreed  quite 
well  with  the  estimate  (Table  2)  but  the  associated  damping  factor  was  found 
to  be  6,7^  from  energy  balance  considerations  and  2,^7’  from  a  polar  plot^S, 
Higher  modes  were  found  at  32.2,  60«0  and  66.7  c/s  whereas  the  estimated  first 
overtone  vms  59.7  o/s.  The  damping  factors  associated  VTith  these  modes  were 
considerably  higher  than  the  assumed  value  of  1/  and  the  amplitudes  at  the 
gyro  station  correspondingly  lower.  For  example,  for  the  fundamental  mode  the 
forced  vibration  amplitude  was  found  to  be  only  15J  of  the  estimate. 


In  Noveniber  1958  revised  estimates  of  the  torsional  mode  frequencies 
etc  were  issued29  as  comparison  of  the  resonance  test  results  with  the  earlier 
estimates  had  shown  unexpected  discrepancies.  These  v/ere  due  to  the  effect  of 
flexibility  in  the  fin  root  fittings  vdiich  made  the  fin  bending  frequencies 
much  lower  than  had  been  anticipated.  Static  tests  on  a  single  fin  attached 
to  a  motor  bay  structure  were  carried  out  to  confirm  this  theory.  The  revised 
estimates  gave  the  fundamental  torsional  mode  frequency  as  56.8  c/s  with 
J4  =  1.09  X  10“°,  podded  fin  modes  at  22,3  and  25.55  c/s  with  gyro  amplitude 
*  —6  —6 

factors  of  0,92  x  10  and  5.60  x  10  respectively  and  a  plain  fin  mode  at 
73.0  c/s  with  an  anplitvide  factor  of  4*64  x  10~6  assuming  critical  damping. 
The  effect  of  these  on  the  modified  roll  stabilisation  system,  i.e,  with  L,C, 
filter  included,  was  examined  but  in  no  case  is  instability  Indicated,  The 
worst  case  occurs  for  the  25.55  c/s  mode  trfiich  gives  a  feedback  vector 
magnitude  of  0.63  at  this  frequency  with  a  phase  lag  of  289°.  It  appears 
hoiiTever  that  the  assumption  of  1^  damping  is  pessimistic.  The  resonance  tests 
did  not  give  resxilts  from  which  the  torsional  mode  damping  covild  be  satis¬ 
factorily  determined^O,  d\is  to  impurity  of  the  modes  excited,  but  the  polar 
plots  suggested  that  it  was  appreciably  more  than  this.  Torsional  resonances 
were  fovind  at  19.7  c/s,  22,8  c/s  26,0  c/s,  45.6  c/s,  54«1  c/s,  64,5  c/s, 

65.4  c/s,  71.2  c/s  \7ith  tanks  enpty  and  19*2  c/s,  21,6  c/s,  48.2  c/s,  55*8  c/s 
and  66,3  c/s  with  tanks  full.  The  modes  with  frequencies  belov;’  30  c/s  are 
predaninantly  fin  modes  and  those  above  50  c/s  have  mode  shapes  similar  to 
the  fundamental  torsional  mode. 


The  results  from  these  resonance  tests  did  not  provide  a  satisfactory 
explanation  of  the  roll  stabilisation  system  noise  troubles  which  had  led 
to  the  insertion  of  the  L.C.  filter.  In  the  hope  of  obtaining  further 
elucidation  direct  measurements  have  been  made  of  the  gyro  demodulator  out¬ 
puts  resulting  ftom  small  amplitude  motor  oscillations  in  a  freely  suspended 
vehicle.  These  measurements  have  been  made  over  a  wide  range  of  frequencies 
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at  1  o/s  intervals  with  tanks  full  and  en^sty*  The  results  obtained  have  not 
yet  been  fully  analysed,  as  they  are  very  conplex,  but  it  appears  that  the 
gyro  outputs  do  not  agree  well  especially  at  hi^er  frequencies  with  what  had 
been  expected.  Possible  reasons  for  this  disagreement  are  the  effects  of 
vibration  on  the  gyroscopes,  local  resonances  in  the  strtxjture  on  which  they 
are  mounted,  play  in  the  gyro  bearings  etc.  It  has  been  found,  for  exaaple, 
that  the  increase  in  roll  rate  anplifier  gain  necessary  to  make  the  system 
unstable  on  the  lavinoher  varies  to  a  marked  degree  fran  one  gyro  to  another. 
Modified  G.¥.5  gyroscopes,  whose  nutation  frequency  had  been  increased  to 
over  100  c/s  by  fitting  a  heavier  rotor,  were  found  to  give  a  larger  gain 
margin,  in  the  absence  of  thmost,  than  the  Series  7  type  used  in  BK.01. 
Subsequent  vehicles  have  therefore  been  fitted  with  modified  G.li7.5  gyroscopes. 

Eight  Black  ICnight  test  vehicles  have  been  fired  to  date  and  in  all  of 
these  the  control  system  has  functioned  satisfactorily  in  flight.  In  every 
case  the  telemetry  records  show  smaDJL  acplitude  pitoli  aiid  yaw  oscillations 
at  about  1*8  o/s  betvreen  75  and  105  sec  after  talco^off.  Those  are 
prestoned  to  be  due  to  propellent  sloshing,  in  view  of  their  frequency  and  time 
of  occurrence,  but  there  is  no  evidence  that  they  are  in  any  way  harmful.  Roll 
errors  have  been  very  small. 

1 0  CONCLUSIONS 


The  analysis  has  shown  that,  in  the  absence  of  propellent  sloshing, 
adequate  control  system  stability  margins  can  be  obtained  with  the  parameters 
selected,  Kence  it  was  concluded  that  the  response  of  the  vehicle  to  guidance 
signals  and  aerodynamic  disturbances  would  be  stable  and  adequately  damped. 
!niis  has  been  confirmed  in  flight.  Due  to  the  low  values  of  danping  factor 
associated  with  the  propellent  sloshing  modes,  however,  pitch  and  yaw  oscil¬ 
lations  of  low  frequency  and  snail  anplitvide  are  likely  to  occur,  anall 
baffles  have  been  fitted  to  part  of  the  oxidant  tank  to  increase  the  damping 
of  these  modes. 

Ground  tests  have  shewn  that  the  higher  frequency  effects  of  struotviral 
flexibility  on  the  control  system  are  not  in  accordance  with  expectation  based 
on  the  simple  theory  used  in  the  analysis.  In  consequence  it  was  found 
necessary  to  modify  the  roll  stabiHsation  system  by  inserting  an  additional 
filter. 


I^p  to  date  eight  Black  Knight  test  vehicles  have  been  fired  and  in  all 
of  them  the  control  system  has  functioned  correctly  in  fli^t.  It  is  therefore 
concluded  that  the  system  is  basically  sound  and  reliable. 

11  FURTHER  DEVllLQEMSM'S 

Development  work  is  well  advanced  on  a  transistorised  version  of  the 
control  system  described  in  this  Report,  j^reclable  saving  of  weight  will 
be  obtained. 

For  stabilising  the  vehicle  attitude  after  the  rocket  motors  have 
stopped  thrusting,  a  bang-bang  jet  control  system  with  dead  space  has  been 
developed.  The  Jet  efflux  will  be  steam  generated  from  catalytic  deoaiposition 
of  hydrogen  peroxide.  This  system  has  not  yet  been  tested  in  flight  but  single 
plane  ground  rig  tests  have  been  successfully  carried  out. 
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LIST  OP  SBfflCiLS 


slug  ft^ 


The  roll  mcment  of  inertia  of  the 
vehicle 


ft 


2 

per  sec 


longitudinal  acceleration  of  the 
vehicle 


slug  ft 


2 


moment  of  inertia  of  the  vehicle  in 
ysvr 


ft  distance  of  the  combustion  chamber 

centre  of  mass  behind  the  pivot  axis 

the  ojrponential  function 


volts 


output  signal  from  the  yaw  command 
vinit 


output  signal  from  yaw  summing  ampli¬ 
fier  variable  gain  potentiometer 


*3 


tt 

It 


f 

^n 

^T1***  ^Tn 
f  p  and  f 


cycles  per  sec 

It 

n 

tt 


output  signal  from  yaw  gyro  demodulator 

output  signal  from  yaw  phase  advance 
amplifier  variable  gain  potentiometer 

frequency 

body  bending  mode  frequencies 

body  torsional  mode  frequencies 

fvindamental  H.T.P,  and  kerosene  slosh¬ 
ing  frequencies 


tt 


frequencies  at  which  the  missile 
transfer  function,  including  propel¬ 
lent  sloshing,  has  zeros  and  poles 

approximate  natiiral  frequencies  of 
closed  loop  heading  control  system 


fy  and  f^ 


G  (p) 


g 


ft  per  sec* 


8nP«  «nL  ‘“n 


n 


freqxiencies  associated  with  the  kero¬ 
sene  sloshing  mode  under  closed 
loop  conditions 

frequencies  associated  with  the  H.T.P. 
sloshing  mode  under  closed  loop 
conditions 

frequency  dependent  part  of  the 
transfer  function  Y^(p) 

acceleration  due  to  gravity 

constants  for  the  nth  bending  mode 
defined  in  section  2 

constant  for  the  nth  torsional  mode 
defined  in  section  8 


i 
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LIST  OF  SYMBOLS  (Cont'd) 


n 


kg 


M 


M 


m 


nig 


N 


N 


N 


n 


PI,  P2 

P 

R 

R 


lb  per  ft 

lb  ft 

lb  ft 

lb  ft  per  radian/sec 
ft 

ft 

slugs 

slugs 

slugs 

slugs 


lb  ft  per  ft  per  sec 


gain  factor,  which  is  independent 
of  freqi.iency,  associated  with 
transfer  function  Y  (p) 

spring  rates  associated  respec¬ 
tively  mth  H.T.P,  and  kerosene 
fundamental  sloshing  inodes 

reaction  couple  at  motor  pivot  due 
to  angular  acceleration  of  ccmbus- 
tion  chamber 

roll  disturbing  torque 

roll  dan5>ing  coefficient 

distance  of  effective  centre  of 
gravity  from  motor  pivot  station 

distances  of  sloshing  masses  m^  and 
mg  behind  effective  centre  of  gravity 

mass  of  empty  vehicle  plxis  ”fi»ed'* 
part  of  both  propellents 

mass  of  one  combustion  diamber 

mass  of  vehicle 

moving  masses  associated  respectively 
v/lth  fundamental  H.T.P,  and  kerosene 
sloshing  modes 

number  of  encirclements  of  the  point 

aerodynamic  yawing  moment  derivative 
%7ith  respect  to  sideslip  velocity 


lb  ft  per  radian  per  sec  aerodynamic  yawing  moment  derivative 

with  respect  to  rate  of  yaw 

suffix  denoting  number  of  bending 
or  torsional  mode 


number  of  poles  of  feedback  transfer 
function  in  right  hand  half  of 
p-plane 

rocket  motors  controlling  vehicle 
in  pitch  plane 

Laplace  transform  complex  variable 
cr  +  jo) 

radius  at  which  motor  thrust  acts 

radius  of  g3nration  of  enpty  vehicle 
plxxs  fixed  part  of  propellents 

rate  of  yaw  (also  =  1) 


ft 

ft 

radians  per  sec 
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T  lb 

Tj  lb 

Tj  sec 

^  S6C 

u  ft  jjer  sec 

u  rad  per  sec 

n 

V  ft  per  sec 

i„(p) 

Y^(p)  =  K^G^(p) 

YgCp)  =  KgCgCp) 

Y^(p)  =  K^G^(p) 

\(P)  =  \\(P) 

Y5(p)  =  k^g^Cp) 

YgCp)  =  %(p) 

Y^(p)  =  K^G^Cp) 

Y(>)  = 

Y^  lb  per  ft  per  sec 

Y^  lb  per  radian  per  sec 

Y1,  Y2 
>  Yg 
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(Cont*d) 

longitudinal  conrponent  of  motor 
thrust  =  cos  5^° 

total  thrust  from  four  motors 

gyro  demodulator  time  constant 

time  (from  launch) 

longitudinal  ccn5)onent  of  missile 
velocity 

2%  X  frequency  of  motor  oscillation 
at  •which  excitation  of  the  nth  bend¬ 
ing  mode  due  to  motor  thrust  is 
cancelled  by  "the  inertia  reaction 
couple 

sideslip  -velocity  in  the  ya-w  plane 

open  loop  transfer  function  (or  feed¬ 
back  transfer  function) 

rocket  motor  servo  mechanism  -transfer 
function 

missile  dynamics  transfer  function 
(yaw  plane) 

gyro  demodulator  transfer  function 

yaw  phase  advance  amplifier  transfer 
function 

yaw  summing  amplifier  -transfer 
fvinctlon 

roll  sunaning  amplifier  transfer 
f\inction 

roll  phase  advance  amplifier  -transfer 
function 

frequency  response  function  corre pond¬ 
ing  to  transfer  function  Y(p) 

aerodynamic  normal  force  derivative 
with  respect  to  sideslip  velocity 

aerodynamic  normal  force  derivati-ve 
•with  respect  to  yate  of  yaw 

rocket  motors  con-fcrolllng  -vehicle 
in  yaw 

lateral  displacements  of  masses  m. 
and  m„  -wi-th  respect  to  vehicle  oen-tre 
line 
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tp 

'I's 
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LIST  OP  SYMBOLS  (ContM) 

number  of  zeros  in  denoninator  of  closed  loop 
transfer  function  in  ri^t  hand  half  of  p-plane 

radians  deflection  of  rocket  motors  used  for  yaw  plane 

control 

radians  differential  rocket  motor  defleotion  giving  rolling 

mcment 

damping  feustcr 

dajiping  factor  associated  with  frequency  f^  etc 


radians 

angle  of  roll  of  vehicle  with  respect  to  roll 

datum 

radians 

an^e  detected  by  roll  gyroscopes 

radians 

vehicle  heading  angle  in  yaw  plane  vrLth  respect 
to  vertical 

radians 

angular  deflection  in  yaw  at  gyro  station  due 
body  flexure 

to 

radians 

angle  detected  by  yavr  gyroscope 

radians  per  sec 

angular  frequency  =  2%  f 

radians  per  sec 

2  f  etc 
n 

NOTATIOT 


Dot  denotes  derivative  with  respect  to  time  e,g.  v  =  tt-  , 

qx 

oo 

Bar  denotes  Laplace  transform  e,g,  =  j  e"*^^  dt. 
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TABLE  1 

Estimated  bending  mode  natural  frequenoiee  and  forced  vibration 

amplitudes  at  gyro  station  for  excitation  at  motor  pivot  for 

riveted  Black  Knight  structure  (G4/2000  Series) 


Missile  condition;  time  of  flight 
on  nominal  trajectory;  bending  mode 

Bending 

mode 

frequency 

c/s 

Angular  amplitude  in  radians  of  forced  vibra¬ 
tion  at  gyro  station  at  resonance  due  to 
Qxcltati'jn  at  motor  pivot  of  (a)  1  lb  lateral 
force  and  (b)  1  lb  ft  couple  assuming  K  of 
critical  damping  =  0.01). 

(a)  Lateral  : 

force  of  1  lb 

(b)  1  lb 

I  ft  couple 

Tanks  full. 

t  =  0  sec 

X  10”*^ 

Fundamental 

(f,) 

7.55 

39.8 

X  10”^ 

5.7 

1st  overtone 

(f2) 

16,5 

11.4 

II 

2,5 

tl 

2nd  " 

(fJ 

25.7 

5.03 

It 

1.27 

tl 

3rd  " 

(V 

32.9 

2, 98 

It 

0.87 

tl 

50^  burnt,  t 

=  72,5  sec 

Fundamental 

(f,) 

9.68 

35.6 

It 

6,23 

It 

1st  overtome 

(fp) 

27.6 

7.56 

It 

2,38 

tl 

2nd  '• 

(fj) 

35.4 

3.75 

tl 

1.53 

tl 

3rd  ” 

47.3 

0.505 

It 

0,271 

It 

68,75?  bijrnt. 

t  =  1CK)  sec 

Fundamental 

(fl) 

12.34 

30.0 

If 

6.49 

tt 

1st  overtone 

(f2) 

33.9 

5.43 

tl 

1.83 

tl 

2nd  ** 

(f3) 

44.3 

1.70 

II 

0,84 

tl 

87*55?  burnt. 

t  =  127  sec 

Fundamental 

(fl) 

18.6 

23.6 

It 

7.42 

It 

1 st  overtone 

(fp) 

41.3 

4.92 

tl 

1.91 

tt 

93*^  bxxrnt. 

t  =  136  seo 

Fundamental 

(f,) 

21,4 

20.9 

tf 

7.84 

It 

1st  overtone 

(fj) 

49.6 

5.07 

11 

2.52 

tl 

All-b\irnt,  t 

=145  seo 

Fundamental 

it,) 

34.9 

11.9 

II 

8.97 

tl 

1 st  overtone 

(f^) 

73.8 

\ 

2.86 

It 

2.77 

tt 
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TABLE  2 

Estimated  bending  mode  natural  frequencies  and  forced  vibration 

amplitudes  of  vehicle  at  gyro  station  for  excitation  at  motor 

pivot  for  Trelded  Black  Knight  struotvire  (G4/3000  Series) 


Missile  condition;  time  of 
flight  on  nominal  trajectory; 
bending  mode 

Bending 

mode 

frequency 

c/s 

Angular  amplitude  In  radians  of  forced  vibration  at  gyro 
station  at  resonance  due  to  excitation  at  motor  pivot  of 
(a)  1  lb  lateral  force  and  (b)  1  lb  ft  couple  assuming 
of  critical  damping  k  0«01) 

(a)  Lateral  force  of  1  lb 

(b)  1  lb  ft  couple 

Tanks  full  (t  =  O) 

Fundamental  (f  ^  ) 

8.17 

35.1  X  10“^ 

5.14  X  10“^ 

Ist  overtone  (fg) 

19.1 

CD 

• 

CD 

3 

2,01  " 

2nd  ''  (f^) 

31.3 

3.27  " 

0.93  " 

3rd  "  (f^) 

42.4 

0,98  " 

0.32  " 

‘burnt  (t  =  72»5  sec) 

Fundamental  (f ^ ) 

10.63 

29.8  " 

5.33  " 

1st  overtone  (f^) 

28,6 

6, 60  '* 

1.75  " 

2nd  •'  (f^) 

39.5 

2.12  " 

0.64  " 

3rd  "  (f^) 

68,4 

0.04  " 

0.03  " 

66.7^  burnt  (t=96.7) 

Fundamental  (f ^ ) 

12,9 

26,2  " 

5.48  •' 

1st  overtone  (f^) 

30.2 

5.36 

1.38  ” 

2nd  "  (f^) 

50,5 

0,60  " 

0,23  " 

83.35^0  burnt  (t=  120.8) 

Fundamental  (f ^ ) 

16.4 

21,7 

5.91  " 

1st  overtone  (fg) 

33.7 

5.36  " 

1.44  " 

2nd  "  (f^) 

80,3 

0,33  " 

0,06  " 

All-burnt  (t=1if5) 

Fundamental  (f ^ ) 

22,1 

9.45 

6,26  " 

1st  overtone  (fg) 

59.7 

1.28  " 

1.83  " 
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FIG.I.  BLACK  KNIGHT  CONFIGURATION 
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FIG.  4.  BLOCK  DIAGRAM  OF  THE  HEADING  CONTROL 

SYSTEM  FOR  THE  YAW  PLANE. 
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FIG.  5. 
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FIG.  5.  DIAGRAM  SHOWING  PROPELLENT 
SLOSHING  MECHANICAL  ANALOGUE. 
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FIG.  6.  ESTIMATED  VARIATION  WITH  TIME  OF 
VEHICLE  MASS,  LATERAL  MOMENT  OF  INERTIA 
AND  CENTRE  OF  GRAVITY  POSITION. 
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RG.7.  ESTIMATED  VARIATION  WITH  TIME  OF  THE  AERODYNAMIC 
DERIVATIVES  AND  FOR  ZERO  INCIDENCE. 
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FIG. 8.  ESTIMATED  VARIATION  WITH  TIME  OF  TOTAL 
MOTOR  THRUST,  VEHICLE  NET  ACCELERATION,  VELOCITY, 
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FIG.  9. 
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FIG.I4. 


TORQUE  MOTOR 


6 


LOAOINQ  SPRING 


ei. 


A.C.  SUPPLY 


FIG.  14.  HYDRAULIC  SERVO 
SCHEMATIC  LAYOUT. 
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FIG.  16. 


FIG.  16.  SERVO  AMPLIFIER  CIRCUIT  DIAGRAM. 
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FIG.  18. 
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FIG.18.  DEMODULATOR  FREQUENCY  RESPONSE  CURVES. 
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FIG.  19. 


/lENTARY  HEADING  CONTROL  SYSTEM  AT  t=80SEC. 
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FIG.  20. 


WG  CONTROL  SYSTEM  AT  t«80 SEC. 
tODES  BUT  WITHOUT  STABILISING  FILTER. 
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FIG. 21.  FREQUENCY  RESPONSE  CURVES  OF  THIRD  ORDER  FILTER. 

TRANSFER  FUNCTION  IS  /,4.o.o.8  aV»: 


GAIN  FACTORS  K,  K3K4=  I  AND  K,  K3  Ks  =  4- 
NOTE.  change:  of  scale:  at  2  0  =  6  DB. 


G.W/R9930. 


P 


U 
y 
o 


U  CQ- 


UJ  KQ 


,  MAGNITUDE  magnitude  I  110.5 


I 

I 


Gw/p,  9931 


REPORT  G  W.  £4-. 

FIG,  23. 
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FIG.  23.  CIRCUIT  DIAGRAM  OF  LATERAL  SUMMING  &  PHASE  ADVANCE  AMPLIFIERS. 
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FIG.  24. 


FIG.  24.  FREQUENCY  RESPONSE  CURVES  OF  THE 
LATERAL  PHASE  ADVANCE  AMPLIFIER. 
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FIG.25. 


FIG.25.  NYQUIST  DIAGRAM  OF  HEADING  CONTROL  SYSTEM 

AT  t=  145  SEC. 
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FIG.  2  6. 


FIG.  26.  NYQUIST  DIAGRAM  OF  HEADING  CONTROL  SYSTEM 

AT  TAKE-OFF  t  =  O. 
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FIG.  27. 


FIG.  27.  PROPOSED  VARIATION  OF  GAIN  DURING  FLIGHT. 
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FIG.29. 
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FIG.  29.  NYQUIST  DIAGRAM  OF  HEADING  CONTROL  SYSTEM 

AT  TAKE-OFF  WITH  REVISED  BODY  FLEXURE  MODES. 

(see  table  2) 
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FIG.30. 
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FIG. 31.  VARIATION  WITH  TIME  OF  FUNDAMENTAL 
PROPELLENT  SLOSHING  FREQUENCIES. 


FIG.  32.  VARIATION  WITH  TIME  OF  FREQUENCIES 
fq.fb.fc  ANDfd  AT  WHICH  THE  MISSILE  TRANSFER 
FUNCTION,  INCLUDING  PROPELLENT  SLOSHING, 

HAS  POLES  AND  ZEROS. 
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FIG.  33. 


FACTORS  ASSOCIATED  WITH  PROPELLENT  SLOSHING 
UNDER  CLOSED  LOOP  CONDITIONS,  NEGLECTING 

COMPONENT  LAGS. 
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FIG.  34. 


FIG.34.  VARIATION  OF  FREQUENCIES  &  DAMPING 
FACTORS  ASSOCIATED  WITH  PROPELLENT  SLOSHING 
UNDER  CLOSED  LOOP  CONDITIONS  INCLUDING 
EFFECT  OF  COMPONENT  LAGS. 
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FIG.  35.  NYQUIST  DIAGRAMS  OF  THE  RIGID  BODY  ROLL  STABILISATION  SYSTEM.  J 
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FIG.  3  6. 


FIG. 36.  NYQUIST  DIAGRAMS  OF  THE  ROLL  STABILISATION  SYSTEM 
INCLUDING  ORIGINAL  ESTIMATE  OF  FUNDAMENTAL  TORSIONAL  MODE. 
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